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Abstract 

We parameterize the phase space density by time dependent clif- 
feomorphic, Poisson preserving transformations on phase space acting 
on a reference density solution. We can look at these as transforma- 
tions which fix time on the extended space of phase space and time. In 
this formulation the Vlasov equation is replaced by a constraint equa- 
tion for the above maps. The new equations are formulated in terms 
of hamiltonian generators of one parameter families of diffeomorphic, 
Poisson preserving maps e.g. generators with respect to time or a 
perturbation parameter. We also show that it is possible to param- 
eterize the space of solutions of the Vlasov equation by composition 
of maps subject to certain compatibility conditions on the generators. 
By using this composition principle we show how to formulate new 
equations for a hybrid fluid kinetic theory. This is done by observ- 
ing that a certain subgroup of the group of phase space maps with 
generators which are linear in momentum correspond to the group 
of diffeomorphic maps parameterizing the continuity equation in fluid 
theory. 

1 Introduction 

The development of collisionless plasma physics in recent years have had some 
interesting breakthroughs. Among these should especially be mentioned the 
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recently found Eulerian action principles for plasma physicsfXarsson 1 ' 2 , Ye 
and Morrison 3 , and Fla 4 ). The intention of this work is to demonstrate 
that some of the substantial problems in developping good models in plasma 
physics comes partially from purely formal problems with how to formulate 
model equations in an efficient , invariant language. A key problem is how to 
discriminate between the incoherent (containing heat fluctuations and reso- 
nant particle interactions) and coherent, fluid part of the plasma fluctuations. 
For this purpose we will use the flexibility of the generator approach to sep- 
arate between generators of the fluid motion to first order in momentum and 
an incoherent part of the generators to higher order in momentum. 

We will be interested in reformulating the continuity equation and the 
Vlasov equation in terms of the action of infinite dimensional transforma- 
tions on space and phase space on the density and phase space density. By 
parameterizing the densities by these transformations with respect to refer- 
ence densities, we find that the continuity equation for densities in fluid and 
kinetic theory are consistent with that the action of the transformations on 
reference fluid velocity and hamiltonian vectorfield plus a timelike generator 
of the transformations are constrained to be equal to the actual fluid veloc- 
ity and hamiltonian vectorfield. We thus replace the continuity equations by 
constraint equations for the time dependent transformations. Moreover, we 
demonstrate in App. A that the continuity fluid equation and the Vlasov 
equation can be looked upon as the defining equations for an infinite di- 
mensional pseudogroup on space and phase space extended with time, but 
constrained to transformations which fixes time. The moment we have re- 
alized this an interesting composition principle pops up. Namely, we have 
not only discovered one new equation, but an infinitely many corresponding 
to different compositions of tranformations which are compatible with the 
above density structures. Our philosophy is that we encode a priori informa- 
tion into the choice of composition transformations. The a priori information 
restrict the class of experiments or processes which the constructed theory is 
intended to describe, and every composition will indeed give new equations. 
In this respect our philosophy is analogous to the theory of measurements 
in quantum mechanics where every measurement is related to a new class 
of operators. In our case, for every measurement or class of experiments 
there correspond a composition principle and corresponding equations. In 
our opinion this principle put the construction of invariant wave equations 
and field theory into a new light (and we believe not only in plasma physics) 
since to every a priori information encoded there correspond new wave oper- 
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ators with corresponding new spectrum and dispersion relations. 

In this work we will try this new principle by making a hybrid fluid-kinetic 
theory for collisionless plasma physics. Thus we restrict our attention to a 
class of experiments where it makes sense to measure fluid density and ve- 
locity. If these quantities are not reasonable to measure, the theory of course 
does not apply. The moment we have decided on this extra bit of a priori 
information, we have a new theory since we can explicitly construct and sep- 
arate out the fluid generator from the incoherent generator containing heat 
fluctuations and resonant particle interactions.^ It must be noted that in the 
ordinary Vlasov equation the heat fluctuations are not explicitly described as 
mode of fluctuation like e.g. in fluid theory. Even in the constant background 
case the heat fluctuations will be hidden in the continous spectrum. With our 
invention of hybrid fluid kinetic theory the heat fluctuations will appear as 
a mode in any background or size of fluctuations due to the new operators 
appearing. 

The plan of the paper is that we will use canonical coordinates to separate 
the fluid generators from the the incoherent kinetic generators. This is done 
to obtain as simple presentation as possible. Then we introduce what we call 
interaction physical coordinates in which the Poisson bracket is noncanonical, 
but still fixed by reference electromagnetic fields . These coordinates are 
essential since otherwise we would have to perturb the brackets also. Then 
we introduce kinetic fluid generators which is related to the above, but which 
can be interpreted as near identity transformations and moreover coincide 
with the fluid generators when we integrate the phase space density over 
momentum. 



1 This notation is however not quite precise since for an experimental situation where it 
is reasonable to define and measure also the stress tensor as a fluid variable, one would try 
to separate out a fluid equation also for this variable by an additional transformation with 
a generator which is second order in momentum. Similar developments and introduction 
of a priori information can in principle be done to any order in momentum and thus put 
tighter and tighter constraints on the class of experiments and measurements which can 
be described by the corresponding theories. The resonant particle effects will for such 
theories be described by higher order generators (in momentum) in interaction with the 
the corresponding fluid variables.. 
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2 Parameterization of the Vlasov phase space 
density 

For a given reference distribution f°, it turns out that the accessible part 
of the space of distributions can be traced by canonical displacements, i.e. 
% = {S,f}, where S is any hamiltonian (e.g. Larsson 1 ' 2 ). These are the 
only allowable displacements in a collisionless plasma. Another way to pa- 
rameterize such displacements are by near identity transformations (see the 
references 3,4). 

/ = exp(£ w )/°, 
C w = {w, ■} . 

As we will see a more general way of expressing the above result is in terms of 
the action of Poisson preserving maps / = ip» f°. (The notation is explained 
in Appendix A and in the text below.) 

In the litterature one has often considered this parameterization of solu- 
tions of the Vlasov equation to be a result of the canonical transformations 
resulting from the underlying particle orbits. A different , may be more nat- 
ural point of view, is to consider the above parameterization as a result of 
that the Vlasov equation is a Lie equation 5,6 having an infinite dimensional 
symmetry group preserving density on extended phase space. In App. A we 
have elaborated on this point of view. 

The Vlasov equation in canonical coordinates has the form 

^ + {f,H} = 0. (1) 

This equation simply expresses conservation of phase space density, u(t) = 
f(z,t)d 6 z, in a Hamiltonian flow. In this report we will basically use canon- 
ical phase space coordinates to derive our hybrid fluid kinetic theory since 
this leads to that the brackets are not perturbed and the Jacobian is unity for 
canonical transformations. This is approach give simple derivations, but leads 
to no loss of generality. We will show how to apply the method in physical 
euclidean coordinates^ also by introducing a fixed Poisson bracket in phys- 
ical coordinates which is not perturbed (see App. A and below). In effect, 

2 Since the fluid and electromagnetic fluctuations naturally are divided into divergent 
and divergcnccfrcc parts it is actually necessary to introduce these concepts on Riemannian 
manifold with respect to a metric even if the laboratory frame is euclidean. The reason 
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if one restricts to canonical coordinates it is possible to use the canonical 
distribution function instead of the density volumeform as the basic entity. 
We demonstrate in App. A that eq. ([!]) is equivalent to that distribution 
function is parameterized by canonical transformations on the phase space 
with respect to a reference solution of the Vlasov equation 

/ = f 1 7^/°of 1 . (2) 

Here ip is a canonical (i.e. Poisson bracket preserving) transformation of the 
phase space P, i.e. ip~ u {g, h} = {ip~ 1 *g,ip~ 1 *h}. The infinitesimal version 
eq. (Q) is expressed by hamiltonian generators. The hamiltonian generator 
with respect to the time parameter is given by 

/,* = /} + ^- 1 *(/°,t ) = M + iT^o, /} (3) 

Here we have assumed that the distribution function and that the canonical 
transformation depend parameterically on t, i.e. f(t) and ip(t). If we in 
addition assume that they depend on one (or several) additional parameter 
e, i.e. f(t,e) and if)(t,e), we can also define a hamiltonian generator, ip e , 
with respect to e as 

f,e = R, /} (4) 

We could think of this additional parameter as a formal perturbation 
parameter which vary say between to 1 corresponding to f°, Hq and f(t), H 
respectively. It could, however, have other interpretations (e.g. describing a 
one parameter symmetry). In the case that the transformation is composed 
of two canonical transformations we have 

•0 = *0 o 0, (5) 
ip* = {jp o ^) _1 * = "0 X * o ^ , 

For composite transformations it is realized that the generators has to 
obey the following rule since they are derived from derivatives with respect 

why is that the above physical division of fluctuations are related to Hodge decomposition 
which transforms in a nontrivial way with respect to diffeomorphisms. We will return to 
this invariant presentation of the fluid and electromagnetic theory in coming papers. 
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to parameters(see App. A) 

A = *Pt + ip *i>u ( 6 ) 

We also know from the assumption that we can interchange the e, t deriva- 
tives for the distribution function that the following Maurer-Cartan relation 
must hold [App. A J 



fite — /jet i (7) 

Vv -i> e ,t +{^ t , i> J + he = o, 

k u = ij~ u kl 

{kuJ} = 0,{A4,/°} = 0. 

Here we obviously can extend our notation by treating the phase space 
coordinates z^ i — 1, ..,6 as parameters and define the hamiltonian genera- 
tors if> i Q (Similarly one can generalize to generators for other parameters like 
the noncanonical guiding center coordinates or oscillation center coordinates. 
See a brief discussion of presentation in other coordinates in App. B). 

fn={AJ} + ^~ U f°n , l'=l,..,6. 

The compatibility condition or the Maurer-Cartan relation then takes the 
form for the seven coordinates (z, t) (here we don't use the e coordinate) 



^i,j-^jn+{^i,^j} + kij = 0, i,j = 1,..,7, (8) 

h- ■ - 7/; _1 *t° h- ■ — —h- 

{%,/} = 0,{*&,/°} = 0. 

3 Notice that the t and e index in e.g.tp t and ip e are not derivatives. Rather these are the 
hamiltonian generators corresponding to the hamiltonian vectorfields ip t and tp e defined 
in App. A. Derivatives will be distinguished from an index by a comma or by explicit 
derivation symbols. 

4 Our definition of hamiltonian generators with respect to different parameter variations 
reflects the fact that our parameterization of the Vlasov equation are Poisson preserving 
maps. Only in the case that the evolution of the background distribution is hamiltonian 
or independent with respect to the parameter in question. Hamiltonian evolution will 
therefore be realized for time and perturbation parameters which we already have observed. 
We also realize that in case the background distribution has an ignorable coordinate- i.e. 
a symmetry, the evolution with respect to this coordinate of the phase space density will 
be hamiltonian. 
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The Hamiltonian generators and the Maurer-Cartan relation can only be 
understood invariantly in the language of forms. We will not describe this 
more general formalism here. 

3 Hybrid fluid-kinetic theory 

In an earlier paper 4 we have elaborated on the formal expansions of the 
distribution function with respect to near identity symplectomorphisms pa- 
rameterized by exponential maps. We purposely did not use the term pull 
back map in that paper, but all the results can be verified by interchanging 
the exponential map with the pull back mapQ. In this work we found the fol- 
lowing equation for the the Hamiltonian generator in the time direction and 
consequently the compatibility condition if we have a additional parameter e 

^ = H-^ U H , (9) 

Let us now use our formalism for composite transformations to try to develop 
interacting equations for fluid and kinetic degrees of freedom. We now assume 
that our distribution are described by 

f = (^°)- 1 7°. (ii) 

Here ip and ip is supposed to contain the fluid degrees of freedom (i.e. mass 
density and momentum) and the 'incoherent' kinetic degrees of freedom re- 
spectively. ip is the map "0 in the reference state. Notice that the order of 
the composition of the fluctuating and the averaged transformations are in 
the opposite order than usually used in passive coordinate transformations 

5 In ref.4 we used partial integration in the variational functional to invert the action of 
the exponential symplectic transformations from one object to another. This trick cannot 
be done with generic pull back maps instead of exponential symplectic maps, but the end 
results are still valid. The reason is that one obtain variational equivalent functionals after 
doing partial integration with respect to the exponential symplectic transformations. We 
define variational equivalent functionals to mean all functionals which gives the same result 
after variation. The trick of partial integration can still be performed after performing 
variations (see App. B). 
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where one want to define hypothetical averaged coordinates, (c.f. Fla 4 where 
we briefly discussed the opposite ordering.) The reason for our choice is our 
goal of separating out the fluid generators from the resonant distribution 
/ and obtain a new Liouville equation for this distribution. With another 
goal in mind other choices could very well be preferrable. To avoid addi- 
tional tranformations due to the perturbations of the generators this is more 
suitable. 

In App. A we have discussed a similar formalism as the above for ideal 
fluid theory and time dependent volume preserving transformations on space- 
time which fixes time (see discussion in App. A). In this case one has to use 
vectorfield generators ij) t and ip e . The parameterisation of the mass density 
and the constraint equation for the fluid generator given reference density 
and velocity are 

u = V>t + ^*u , (12) 
p,t = -V • (up) , 

P , t = -v-(<M 

5 P = -V-(cWv), (13) 

*t = d -^o^-\^= d -±o^~\ (14) 

5^ = Si/joiI?- 1 . (15) 

The parameterized velocity field u give a constraint equation for the above 
diffeomorphisms when we give the velocity field u = u , e.g. as in our case 
from the momentum equation. Here the invariant object is not mass density, 
but the density form u = pd 3 xdt which the pull back map acts properly on. 
The volume density preserving pull back maps can also be presented as an 
action directly on the density by taking into account the Jacobian of the 
mapping. We introduce the symbol ip» for this action which also can be 
parameterized by a volume density preserving near identity transformation 
in the following way (see App. A) 

P (X, t) = If, • P °(X, t) EE p (^-\ X ,t),t)J . (16) 

In the case of near identity transformations one can further express the action 
in terms of the near identity generator w as 

</>.(•) ee exp(-V-(w.))(-), (17) 
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V>* = exp(ad(w)), (18) 
if) = exp(w). 

Here ad(w) = [•, w] in terms of the standard bracket for vectorfields(see 
definition in App. A). The proof of these expressions follows from using 
■0 = exp(w) in the above definitions of the actions of diffeomorphisms on 
densities and vectorfields and compare it with the action of the above op- 
erators. Notice, that the fundamental operators in the exponential is Lie 
derivatives of the corresponding object with respect to the near identity gen- 
erator vectorfield. Further, we can express the vectorfield generators in terms 
of near identity generators as (see the definition in App. A) 



d dw 
*l>t = (^*-Id)—=iexp(ad(w)) — , (19) 

d dw 
tPe = (^*-Id)— = iexp(ad(w)) — 

Sip = iexp(ad(w))5w, (20) 

. . exp(x) — 1 , . 

iexp(x) = FV ' . (21) 

Ob 

If we use the identities established in App. A to transform from the above 
type of expressions on phase space and Hamiltonian vectorfields to Hamilto- 
nian near identity generators, w, we obtain that 

(22) 



i, 


= exp(X lu ) , 




= exp(-£(X w )) , 








= iexp(C w )— , 




tr \ dw 


= iexp(C w )—, 




= iexp(C w )Sw 




= r 

>~"W • 



(23) 



Here, the Lie derivative C(X W ) = —C w ; act as the operator X w = —C w on 
functions, but act of course differently on other objects. The above inter- 
pretation corresponds exactly to the point of view we proposed at an earlier 
stage in Fla 4 . 
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The above formal expansions in terms of near identity generators can 
be thought of to correspond to Larssons perturbation expansion when an 
ordering is given to operators (see App.C and Larsson 1 ' 2 . However, because 
of our group composition concept, we have considerable freedom when it 
comes to modelling of specific physical processes. At all steps in our theory 
it will be possible to specialize to near identity generators, but we will not 
stress this below. 

We can also think of the above maps as a family of diffeomorphisms on 
space parameterized by time. Here J is the Jacobian of the transformation. 
The negative sign in the near identity transformation is used to obtain an ad- 
equate sign in the continuity equation. We used the same reason for positive 
sign in the parametrization of the phase space density which corresponds to 
negative sign in the corresponding Hamiltonian vectorfield. The composition 
of the above maps is similar as for symplectic pull back maps. The compat- 
ibility condition with respect to an additional parameter e, is given by (see 
App. A) 



A,t-^t,e-[^a^t]+^t = 0, (24) 



V-(k eiP ) = 0,V-(k>°)=0, 

Kt = tXt- 

Since the velocity field can be parameterized as in eq. (|12|) , we obtain 
that the compatibility relation can also be written as 

^ M -u, e -[^ e ,u] =0mod(kJ . (25) 

It turns out that it is possible to lift the fluid generating maps to kinetic 
theory and consider them as a subgroup of kinetic, canonical transformations 
through the definition of the following canonical generators 6 

V'tfo*) = <MM) +^2t( x ^) = P-u + V 2 t , ( 26 ) 

^2t = - m y-^ A " G ' ( 27 ) 



6 These symplectic transformations cannot be defined through near identity transfor- 
mations since the resonant particle distribution will be following fluid orbits. This is done 
for mathematical convenience of the separation procedure. Later on when we linearize the 
hybrid fluid kinetic equations, we will relate these symplectic transformations to the fluid 
generators with respect to a reference state. 
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1p e (z,t) = 1p le (z,t) +V 2e (x,t) = P- ll) e + i>2e , 

i) u = p • U , ^ u = p • V e • (28) 

Note that we have assumed that the generator ip 2a ,a — t,e depend only on 
x, t. With the above definitions one easily convince oneself that the compat- 
ibility conditions for the barred symplectic transformation, ip la , a — t, e is 
consistent with the fluid compatibility condition in inner product with p. 

The explicit form of the above kinetic fluid generating maps can be de- 
scribed in terms of cotangent lift ^f(t)(t) = T*<f>(t) in composition with the 
fibertranslation by an exact form described below.. In Abraham and Mars- 
den it is proven that the cotangent lift is a symplectic map which preserves 
the canonical oneform 9 on phase space T*M (= p-dx in euclidean canon- 
ical coordinates). For a point a q G T*M = P and v eT^-i^M one find 
that T*(f>(t)a q (v) —a q {T<f){t)v). Define the coordinate functions on T*M by 
Tip : a q >— > -np(a q ) = z = (x, p). The action on coordinate functions is there- 
fore (T*(f>(t))*7rp which for a euclidean metric means that the action of the 
cotangent lift is = Z = T*(f>(t)(z) = (</>(i) _1 (x), (V0(t)(x)) • p). A 

simple way to find invariant properties of the above map is to use that the 
canonical oneform is preserved under cotangent lift, i.e. (here we treat p iy x % 
as coordinate functions) 

9 = %(t)- u e, 

pjdx j = Pjd^(t)- l *x j , 
P 3 = *(t)- u Pj = ((V<l>(t)(x)-i>) r 

Another map which will interest us is fibertranslations by oneforms A^\t) 
on M , defined by ^ A u (t)9 = 9 - n* M A^(t) where vr M : T*M -> M. 
In this case the map is only preserving the symplectic tensor uj = —d9 
if the oneform A^(t) is exact since ip A u (t)uj = uj + ir* M dA.^ (t) . In eu- 
clidean canonical coordinates the fibertranslations correspond to the map 
ip A : (x, p) — > (x, p — A(x,i)) which up to a numerical factor is the transfor- 
mation to euclidean physical coordinates. The fibertranslation give rise to a 
generatorvectorfield on phase space which can be described as Xf- = — ^ in 
euclidean coordinates, but has an invariant description given below (put the 
fourth component of the four oneform to zero and specialize to Euclidean 
coordinates in the invariant description). 

If we extend the canonical oneform to a oneform 6 on extended phase 
space T*X , X = M x R , 6 \t*m= 9 , the extended map ip which fixes time 
will generate a component in the time direction. 
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Theorem 1 ^ 6 = + dS + J t dt = Q + dS + ^ t dt, 

at ri 

ip~ u Q = O - dip t A dt ,Q = -dQ, 



§ = Jt-A- 



Proof. The proof of the above lemma is simply by observing that i^e_ +x ^ ^(ip 

since the extended oneform by definition has no time COm- 
df 

ponents. Similarly, the above relation between the generators and the gauge- 
fields is a consistency requirement coming from £(J^)('?/>~ 1 *fi) = — C(X^ t )(ip~ u n) 

-£{X A )n = -d(#t) 

= -£(f)# t a dt = -c(f t )d(j t -§)Adt.m 

In coordinates J t = P - ^~- in the case of cotangent lift we find that 

Jt = P-^ = % = P ; 0t , 4»t = (0* - Id)M 

The above extension of the canonical one and twoforms to extended phase 
space, give an opportunity to consider the action of noncanonical transfor- 
mations on these forms. Therefore we can e.g. consider transformations 
generated by oneforms on extended space X. 

e A = ^ u e = e-n x A^, 

n A = ^ 1 *n = tt-ir x dA {1) . 

The generating vectorfield corresponding to the above generalized fibertrans- 
lation is given by Xf = -J A • ^ x iX (l) ) , Y {1) = i±(dA^) where J A = ip A J 

dt 

and J is the Poisson tensors corresponding to uj a = £l A \j>* m and uj = Vt \t*m ■ 
In the special case of translations by an exact form = df (here / is a 
function on space/time not phase space,we find that the generating vector- 
field is hamiltonian and the transformation trivially symplectic, i.e. Xf = — 
J-(df )t ). In fact, a large class of transformations which is not necessarily pre- 
serving the the symplectic tensor, but is still divergencefree, is generated as 



7 We can consider our results for J t to be an extension to infinite dimensional pseu- 
dogroups of the momentum map for finitely generated groups and Banach Lie groups (see 
e.g. Abraham and Marsden). Earlier results has been directed towards reduction while we 
concentrate our efforts towards the composition principle for in principle infinitely gener- 
ated groups. Similar results with respect to reduced actions also holds for pseudogroups, 
but we will not study it here. 
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X t a = 3 a ■ (A a ), 3 a = ^ aif 3 ,Q, a = ip a U ^ and A a is some oneform given on 
T*M. 

It is possible to extend the above theorem to a more general setting 
based on noncanonical twoforms VL n = V^ 1 *^ an d oneforms 6 n = V^ 1 *© by 
introducing the convective derivative = + X™ , X™ = — /d) . 

Theorem 2 ^ n )~ 1 *e n = 6" + dS™ + J"cft, 

(V> n ) -1 *^ n = n n - d{i) n T - ^) a dt ,n n = -dQ n , 
ij.((r)- u e n ) = -i x ^((r)- u e n ) = J n T , 

dr ^ T 

dS n _ jn _ i n 
dT ~ T Vt ' i 

j n = 4 n o4oK\X+n = (tl> m -Id)£,S» = rl>; u S. 

We are now in a position to explicitly describe the above fluid kinetic 
generators in canonical, euclidean coordinates as the composition of the fol- 
lowing cotangent lift and translation by an exact form df , i.e. %p = o ip d j 
with generator xjj t = * t + # _1 */,t = 4>t • P+ft , ft = ^~ u f,t • Here we must 
choose t = u = ij> t + V'* u o which for u = Vot defines the unperturbed 
4> such that <p — ij) ° Vv The gaugefunction has to be chosen as = ^ 2 t to 
agree with the above. Since we have assumed uo to be unperturbed, we find 

that fa = V e • P+/e , fe = V- U f,t. 

We now have the following theorem for the action of fluid kinetic maps in 
canonical coordinates defined as a composition between a cotangent lift and 
translation by an exact form 

Theorem 3 Let f — ijj U f such that p = J fd 3 p and p = J fd 3 p . Then 
one has that 

P = 4> • Po = il> • Po , Po = V>o • Po ■ 

Proof. The proof follows from that p = J f ^p*d 3 p = 

J O -1 ) / / d3 P =J (0 _1 )</> _1 *A) = <t>*~Po = *!> • *l>o • Po = 1>*Po- ■ _ 

With the above relation between the Hamiltonian generator and ip t and 

macroscopic fields we obtain the reduced Hamiltonian which can be identified 
with the Hamiltonian in the fluid reference frame. We also find a reduced 
Vlasov equation for the Active phase space density / 

H = H - i) t = (p - muf/2m (29) 
P = P--A, 
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dt 1 U '" J 

H = i?H . (30) 

The motivation for the special choice of ip 2 t{ x J) * s to obtain a reduced Hamil- 
tonian of the above form corresponding to the fluid reference frame. 

We will now discuss the physical consequences of the compatibility con- 
ditions for the distribution function and mass density. The fluid density is 
given by p = J fd 3 p such that the relation between the fluid and kinetic 
generators becomes 

P,t = J(& t J}+V U %)d i p = -V-(up) + jT U ^d 3 p, (31) 

P,e = J /} + ^*%)d 3 P = ~ V • (iM + / T U8 -f/ V . 

Here we have used the fact that / {0(x, t), f(z, t)}d 3 p =g • / ^^d 3 p =0 
with suitable decay properties on / for large momentas. We immediately con- 
clude by comparing with the above fluid theory that the following additional 
restrictions on the hypothetical density / has to be imposed if the barred 
symplectic generators should correspond to a fluid subgroup with respect to 
its moments 

T U %*T> = 0, (32) 
T u % d 3 P = 0. 

It is now obvious that we could add any additional term <f> t (x,t) to the 



generator ip t without changing the form of eq.(31). The trick of using the 
above form of barred generators is however that it is possible to obey eq. 
([32] ) easily. 

We find the following theorem 

Theorem 4 With the earlier definitions one obtains f ip U ^d 3 p =0. More- 
over, the definition of the fluid density from the composed fluid density is 
compatible with fluid theory if we select f such that J ip %d 3 p =0. If we 

select the reduced distribution function such that J if; U ^d 3 p | t= o =0 and 
p(x, t) \t=o= P°(x) , it will also be an identity at any other time. 
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Proof. We need the following lemma 

Lemma 5 For symplectomorphisms, ip, with generators which are not more 
than linear in momentum the following applies for a phase space density g 
with suitable decay properties for large momentas and J *§f = 

g = J tp U gd 3 p = $ • g , (33) 

9 = J gd 3 Y> ■ 

Proof. In the case of near identity transformations the lemma can be proved 
by using near identity generators which are linear in momentum. Since it 
follows that the deviation from identity is also linear in momentum, one find 
the above relation after integration over momentum. If one does not assume 
near identity mappings the proof is a little more involved. The infinitesimal 
version of the above equation is since the generators are supposed to linear 
in momentum and with suitable decay properties of g in momentum such 
that the following holds 

9(t),t = j {^T l \t)g}d*p+ J ^f/v 

= -V • (J uT 1 \t)gd*p)+ J r U %d*p, 
g(t), t = -V-(u$(t)). 

Therefore if j> §f = 0, we can parameterize g(t) = ip(t)»g. With 
this parameterization it also follows that '»/>»§f = from which we find that 
H = since -0 is invertible and takes zero to zero. ■ 

Remark 1 The above parameterization is performed with respect to a ref- 
erence g such that || = 0. We can change this reference by doing an ac- 
tive transformation (or alternatively a passive coordinate transformation) 
g = (00 X ) • 9o su °h that ^ + V ■ (u (?o) = 0. Now it is possible to de- 
fine transformations = ?/> o O 1 which is more suitable for near identity 
transformations with respect to reference state. This is what we will eventu- 
ally do in section 3.2. 
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Using the above lemma we immediately find that 

^•p,a , a = t,e , (34) 

Now we use the reduced Vlasov equation and find 




By partial integration in momentum we find that 

/ r*f i 3 P = -V ■ (/ S (35) 

since by definition f (p — mu)/d 3 p =p(u — u) = 0. (I.e., we naturally will 
restrict the parametrization of the velocity to be equal to the velocity moment 
of the distribution. In fact this restriction is the defining equation for the 
volume density preserving transformation.) 

The compatibility condition for the mass density parameterized by the 
above composition of symplectic transformations acting on phase space den- 
sity leads by eq.'s ([H], |32|) to the constraint 

d 2 p d 2 p 
dtde dedt 

d . f—-i*df ,o . 

=* at J* T/^ 

Therefore we conclude that J ijj ff^ 3 P — V' • / % d 3 p is constant with 
respect to time. Consequently if the e parameterization is such that the 

/ i/j §^ 3 P I £=o = 0' ^ "will continue to be so at all times. We assume that 
both the transformations ip -1 and ip exists. Since the action of the sym- 
plectomorphisms is such that it preserves volume density forms, it is realized 
that the zero density must be transported to zero density by the action of all 
symplectomorphisms. It follows that J ^d 3 p =|£ = if it is initially chosen 
in such a way. Moreover, it is implied that p = p° is given by the reference 
density independent of e even if / is depending on e. By a similar argument 



— -1* 



1p "f,ad 3 p 



(36) 

= 0. 
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we deduce that = 0, and consequently corresponds to a spatial reference 
density given in the fluid frame of reference. 
■ 

The interpretation of the above result is that to obtain a composition 
symplectomorphism where the volume density preserving transformation is 
described as a subgroup of the group of all symplectomorphisms lead to that 
density perturbations are parameterized in phase space by the fluid subgroup 
consisting of the barred symplectomorphisms. 

Let us now study the continuity and momentum equation more explicitly. 
The Vlasov equation can be written 

% + {f,-$t + B) = 

H = -$ t + H. 

The zeroth order moment integrated over momentum space now gives the 
momentum equation with no contribution from the {/, H} term. The mo- 
mentum equation can now be found as 

9 J pfd 3 p+ J p(-{ft, /} + {H, f})d 3 p = , 



dt 



(37) 



^(pu)+V-(puu+P)-pf L = 



f L = — (E + u x B) 



m 



P 



f—(p p -mu)(p p -mu)fd 3 p. (38) 
J m 



Here f& and P are the Lorentz force and the stresstensor of the fluid and the 
physical momentum is related to the canonical momentum by p p = p— |A. 
The reduced Vlasov equation is given in eq. ( p9|) and form together with the 
continuity equation and the above momentum equation a new set of equations 
for collisionless plasma physics, the hybrid fluid-kinetic theory. 



3.1 Parameterization of the hybrid fluid-kinetic theory 

Let us briefly discuss the parameterization of the continuity equation and 
the reduced Vlasov equation. We recall that fluid generating vector is re- 
lated to the parameterized velocity by u = i/; t + V , * u o- We purposely have 



17 



been operating with the parameterized velocity u generated by the diffeo- 
morphisms different from the velocity vector derived from the moment of the 
distribution. The reason why is that they are not a priori equal. Indeed, 
the equality of these two quantities is the constraint equation which together 
with the above fluid momentum and reduced Vlasov equations determines 
diffeomorphism 



u = u . (39) 

This equation replaces the continuity equation since the density can imme- 
diately be mapped the moment we have e.g. the near identity representation 
of the diffeomorphism. The reduced Vlasov equation can be parameterized 
in a similar way as the Vlasov equation by 



4> t = H-^~ U H , (40) 
H = H + ^° t . 

If we parameterize our diffeomorphism by the perturbation parameter e, we 
use the compatiblity equation to obtain the the determining equation for the 
generating function ip e 

^-{#,^}-tf e = 0. (41) 



Eq.'s (|37|-[41|) now form a complete set of equations as an alternative to the 
continuity and reduced Vlasov equation together with the momentum equa- 
tion. 

We now have to discuss the momentum equation more carefully to ob- 
tain an invariant description of the parameterization. To obtain such an 
invariant description it is useful to describe the velocity field as a oneform 
and the stresstensor as a symmetric, covariant twotensor P^ 2 -*. It is now 
possible to parametrize the velocity field by using Hodge decomposition with 
respect to the threedimensional metric g to split it in a rotational and di- 
vergent part, u*- 1 -* = —dr] + * s dA^ = u£ + u£ . Here we have used the 
Hodge decomposition with respect to the transformed metric (see below) in 
agreement with the philosophy that the decomposition should be given with 
respect to the standard metric when pulled back to the reference level, i.e. 
uW = i^(t)~ u u^ , uW = — dfj + * go dA.lP. It can be useful to further de- 
scribe the rotational gaugefield by Pfaff decomposition as = adf3 + , 
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where the gaugefield does not depend on the gaugepotential 7. We find this 
parameterization more convenient than the standard Clebsch decomposition 
which does not separate into rotational and divergent parts and does not 
use the metric structure. To split the equation with respect to the above 
structure we write the momentum equation as 

(| + £(u))u« = -I^pW + d^uWj + fW, (42) 



P( 2 ) = J(t)p( 2 ', PfE^^)- 1 *^ (43) 
p( 2 ) = VW _1 *P (2) = J2 P ? ®d(ii){t)- u X] ) (44) 

3 

g = *l>(t)- u g (45) 

Pit) = m- u p (t), 

This equation can also be written in an invariant way as an equation for the 
momentum M =lu ® and/or written in terms of covariant derivatives. 
In addition we can embed this it in extended space by taking the wedge 
product with dt as we did for the density form, but we prefer to defer this 
formulation to another paper. Here we have used that D s (-) = -j^D So (J(t)- ) 
to obtain an invariant form of the equations suitable for transformations by 
timedependent diffeomorphisms. We need the following natural properties 
of the diffeomorphism action 

(| + £(u))o^(t)- 1 * = ^o^ + ^uo)), (46) 



*)) 


o^(t)- 1 * 




*g 




= *j>(t)- u * so , 


d 




= ij>(t)- u od, d< 






= #)- h oi fi , 




u 


= ip(t)*u, u = % 




u 


= V*^o , u = u+ 




£ g (p (2) ) 





-1* / —1* 

■ " ' U I 

-1 * 

o % 



o d, 

1* „/.— 1* 

"«0 5 



d d 
>(2)) . (47) 
Here tp is as before the extension of the timedependent map ip(t) to maps in 
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four space which fixes time. We now observe that our momentum equation 
can be pulled back to the reference level where in fact the Lie derivative 
and the fourdimensional interior multiplication will be linear operators with 
respect to the background velocity field and all quantities tranform in a 
natural way. 



(| + £(u ))u« = -^D^+diU^+f^, (48) 

fl 1} = -i o F( 2 ) = -^(E« + * g0 (u«AB«)), 
F {2) = ^-upm ; E d) = ,/,( t )-i* E M B « = tp(ty l *B^, 

B (l) = +gB (2) = HdA (D ; B (D = * goB (2) ? 

In addition after the pullback all operators and physical quantities according 
to the above will be specified with respect to the background metric go- Our 
philosophy is that all equations and physical fields should be represented 
such that they can be pulled back to the reference level. The pulled back 
equation can now be compared with the equation for the reference solution 
and an equation for deformations from the reference fields can be formulated 
in a strikingly simple way in which many of the terms are linear. This make 
our theory especially attractive from a perturbation theory and complexity 
point of view, but we believe this also has implications for the interpreta- 
tion, predictions and formulation of measurements for physical fields. E.g., 
according to us a linearized equation and fields pulled back to the reference 
level in no way is linear at the original level which even obtain an induced, 
nonlinear metric structure. 

We now introduce the deviations at the pullback level of the physical 
fields from the background fields. Since the equation for the reference fields 
also obey eq]4^ taken at the reference metric go, we can subtract the back- 
ground equation from eq.48. We then obtain an equation for the fluctuating 
quantities suitable for perturbation theory 

(|+£(uo))u« = -l^P^ + i^uW-^ui^ + fl 1 ), (49) 

u« = u« + u«, P( 2 ) = P( 2) + P( 2 ), F^ = F^ 2) + F^\ 
fj, = iuoF^ \ e.t.c. 

Notice that even if we were only interested in linearization and linear quan- 
tities at the original level, the distinction between the pullbacked equations 
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and the equations at the original level will still be essential since the fluctua- 
tions in the pullbacked metric and the pullback map itself will not affect first 
order quantities, but they will affect background quantities to linear order 
presented at the original level. 

Let us complete our discussion of the momentum equation by showing how 
elegant it separates into equations for the acoustic and rotational potentials. 
We split the right hand side of eq. ( j46lj48|) into rotational and divergent parts 
by using Hodge theorem to define the potentials 

-dK + * 8 dR m = f} x) - -DgP®, (50) 

-dk + * so dK {1) ee f| 1} - —D g0 PW. 

Pa 

If we now take the exterior derivative of eq. (f46|) , we obtain the vorticity 
equation 

(^- + £(u))du« = d* s dR w , (51) 
U W = * g dAj 1) = * s (da A dp). 

In fact, up to a potential we can even give the equation for the rotational part 
of the velocity oneform(relative to the fluctuation metric) or alternatively the 
vorticity twoform, 7r^ as 

(| + £(u))u« = vffiW, (52) 
(|- + £(u))7r (2) = d* e dR w , (53) 

7T< 2 > EE du«. (54) 

Moreover, if we subtract this equation from the original momentum equation, 
we otain an equation for the acoustic potential 

(^ + £(ft))i ? = «-i.W 1) - (55) 

All this equations transform in a natural way with respect to diffeomor- 
phisms, i.e. we can obviously formulate pullbacked equations and equations 
for fluctuations with respect to corresponding background. 
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Let us now describe our diffeomorphisms more explicitly and relate them 
to the above potentials. We do this by studying the constraint equation for 
velocity vectorfield reformulated as a constraint on the momentum lifted to 
a oneform; pu^ = pu^\ The compatibility equation for the perturbed vec- 
torfields now become a relation between oneforms(lifted by the fixed metric) 

(p^), t -(puW), 6 + * g0 d(* S0 (p^ A u«)) = 0, 
which can be pullbacked to 

(Po^ { e),t -(p u (1) ),e + * g0 d(* go (Po^e ) A * (1) )) = 0. 
If we now take the divergence of the first equation, we obtain that 

(*g d* So (pip { i ) )),t= -P, e t = (* So d* S o (Pu (1) )), e . 

This is in fact the perturbed, parameterized continuity equation. 

We can also study the compatibility equation with respect to the oneform 
Ju« = Ju«. In this case we obtain 

(j4 1) ), t -(Ju( 1 )), 6 + * go d(* go (M 1) Au«)) = 0, (56) 

(*s d* So (Jtp { e)):t = -J, e t= (* go d* go (Ju (1) )),e- 

Here we have used the continuity equations for the Jacobian J, t +V- {Jif> t ) = 
0, J, e +V ■ (Jip e ) = lifted up to forms. In fact it is possible to find 
material coordinates such that the Jacobian and the density is equal up to a 
constant density, i.e. p = J* Pq. The parameterization of the velocity oneform 
if written with respect to reference metric is given by 

u (1) = -dr] + jg ° go 1 ° *go^ A r 1) - 

This means that Ju^ will only contribute to the above divergence term in 
the continuity equation for the Jacobian through the acoustic potential. If 
we want a parameterization where the rotational part does not contribute in 
the divergence term of the mass continuity equation, we have to use Jm, 
. Anyhow, we pull back the velocity oneform and relate it to the generator 
as a oneform at the reference level and find 

V>f ) = -dfj + * go dA« - u 1} = -dfj + * go dA«, (57) 
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Here we have to choose which generator is first and last corresponding to the 
composition of e.g. i/> = tp r o tp d or vica versa. 

The above decomposition into divergent and rotational generators shows 
us the need to close our system of equations by representing the generators ac- 
cording to the type chosen. Obviously, the divergent diffeomorphisms should 
only have one parameter related to the density structure while the rotational 
diffeomorphisms should have only two degrees of freedom related to the vor- 
ticity structure. We notice that if e.g. if) = exp(eif> 1 ), the generator i/) e = xf) 1 
and similar relations can be worked out for other perturbation parameter 

relations. Therefore the key is to study the decomposition = -0 e + i/? e 
represented in terms the density and vorticity structure respectively. We 
define the density structures u = pdVo, u), e = p, e dVo, u e = p e dVo = if)*u), e . 
The pullbacked perturbed density structure is therefore explicitly related to 
the density perturbation as p t = if)*(^f). Above we have decomposed the 
velocity oneform into a divergent and rotational part with respect to an in- 
variant volume element. Here we find it more natural to consider the dual 
decomposition of pip e , pift t and pu (corresponding to uj ® u e.t.c). In terms 
of the perturbed, pullbacked density structure we find that we can represent 
PoV'e = V'*(pV'e) as ^ t 'go(PoV'e) = ~Pe- Therefore we have that modulo a 
rotational part we can define p^ t = — dV g ^p e and the rotational part can 

be represented as p tj) € = * go dA t . In terms of this parameterization the 
divergent part of the compatibility equation which we found to be equivalent 
with the perturbed continuity equation can be formulated at the pullbacked 
level as 

(^ + £( 3 )(v ))p £ + ^ g0 (v £ p ) = 0, (58) 

v e = V^Vje or = V*v (1) , e • 

Equation's (|58|) and (??) (can be given in perturbed form, but that deserves 
a separate study) contains the description of the acoustic mode and the 
interaction with the rotational mode and the kinetic fluctuations. In the 
linearized case this correspond to a longitudinal wave equation with kinetic 
and rotational effects in any background fluid and kinetic state. What has to 
be done is a detailed study of the deformation properties of the stresstensor 
discussed below with respect to the divergent and rotational diffeomorphisms 
and the incoherent kinetic transformations. 

Let us now discuss the rotational mode in more detail. We define p = 
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tp d ° Pr , Pr = i> r ° Po an d + V • (u r p r ) = , u = $ + ?/>*u r , u r = 
tj> r t + , 0I u o = u c + "0I u o > u c — ^*Uc- We now find the following equiva- 
lent compatibility equations for the rotational generators with respect to the 
pseudogroup connected to the density p r 

V'g,t-u r , e +[u r , Vg] = 0, (59) 
equivalent to xj) r e , t — u c , e +[u c , ipt r e ] = , 

or (p r jpl (1) ), t -(p r iiW), e + * g0 rf(* g0 (p r ^e (1) A uf')) = , 

or^.-uW. + ^^^'Aui 1 ')) = 0. 

Here g c = i/' r *~ 1 go is the fluctuation metric with respect to the rotational 
part of the diffeomorphism in contrast to the total fluctuation metric used 
earlier. 8 The above equivalence is due to that (t/^Uo), e = [V>* u oj V'el which 
is valid in general for any background vectorfield. The last equality is valid 
assuming by definition that div gc (u^) = div gc (i/^ r ^) = 0. In App. C we in- 
troduce a description of the rotational vectorfield on a family of level surfaces 
in threespace given by (5 defined as 

X^ (1) = * s J(ad(3), (60) 
X^ 1} = * s J(aJp), 

= * s AaJ{3), (61) 

g c - 1 (* gc rf(a^))e.t.c. (62) 

We define a new roational bracket structure with respect to vorticity 
situated at the foliations in threespace of (3 for given metric g and invariant 
volumeelement dV by 

{a, ftp dV = * s J(ad(3) A * gc df = <df, > dV, (63) 
{aj}p = Xi(f)^* Sc (* Sc (daAdp)A* s Jf). (64) 

In App. C we prove that for purely rotational vectorfields one have that 

PCx gJ = x! {QiiQ2b , (65) 

8 All the quantities and equations introduced here related to the rotational part can be 
transformed by acting with ijj d to obtain the actual measured quantities. The reason why 
we do this pullback by xp d is that many of the relations we present will be very complicated 
without it. But of course all the relations and quantities can be transported to the real 
measured quantites by the inverse action. 



c 

ir(l) 

^ = 
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X£(/(/?)) = {<*,/(/?)} = 0. 



(66) 



This show that one can think about the rotational bracket as a noncanonical 
Poisson bracket with functions of type /(/?) as Casimirs. 

From the above splitting of the constraint equation for the velocity field, 
we find for the rotational part u c = xh\ + iplu 0c = (ip^ o <fi Qif — Id)jL, u 0c = 
(0o* — ^)§t ' Here we will have one representation for each component 
in the constraint equation. A different representation which concentrate on 
deformations of the rotational bracket structure of the reference state is u c = 
= Xg° + Xj = u° + ^* o (fa - Id)u 0c + u c ,f3° = j,- u /3 
(faofa-Id)^ ee (fa-Id)^r ee fa c = X|° c such that u°= Xf = ^ Tc + 

$*u 0c . Here ^ ee f +u c , G c ee (fa-Id)§- t = fa = Xl-fao(fa-Id)uoc . We 
parameterize the reference rotational velocity as uo c = Xa° and find that 
X p _ -_ u = , i.e. a = ib T + xb a mod(f(/3 )) . The perturba- 

a~ip Tc -ip a c ■ 

tional aspects of this constraint equation with respect to a parameter e can 
be explored by defining (fa o 0, - Id) A = _ / d )-| X £, ^ ee 

^ + </> g , g ee (fa — Id) -fig . To avoid perturbations in the bracket struc- 
ture we should pull back the constraint equation for the rotational potential 
and obtain a = ip Tc + a , a = xp a, x]) T = xj) ip Tc . Here we refer the 
potentials to the bracketstructure derived from the background foliation j3 
and corresponding vectorfields X^°. We therefore find that analogous to the 
theory we have developped before for density equations 

fa,t = [0 e ,u c ] + u c , e , (67) 

V'g ee $~ *V> e 

Equation's ( |6"7D and ( |52] ) (can also be perturbed with respect to e ) give now a 
description of both the vorticity structure and the related diffeomorphism in 
interaction with kinetic fluctuations and the longitudinal fluctuations. The 
above equation for the vorticity (|52]) with zero right hand side for the purely 
rotational case is the direct generalization of the potential description of 
the vorticity equation on fixed two dimensional surfaces to convected level 
surfaces foliating threespace defined by (3° and the rotational potential a . In 
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this case u c = tr = Xf and (f> = Id since the rotational equation becomes 
the defining equation for a Lie pseudogroup related to the rotational structure 

(| + £(u))u« = 0, (68) 
(| + /;(u)y 2 ) = 0. 



The rotational diffeomorhisms are still described by eq. |67], but with triv- 



ial (f> = Id. In this case the vorticity equation is analogous to the Vlasov 
equation with a playing the role of the hamiltonian given by ^tt^ = 
* g0 dx^ 0< ' 1 ' ) = d * g0 d(adj3 ). If we invert this operator we find an expres- 
sion for the rotational potential a as a functional of the vorticity analogous 
to noncanonical hamiltonians on constrained surfaces in classical mechanics. 
This also explains why this type of models pops up in many applications 
in fluid and plasma physics. Notice that our composition principle for the 
velocity induces in principle a finite or infinite series of diffeomorphisms and 
divergent and rotational potentials and a corresponding series of longitudi- 
nal and rotational compatibility equations and bracket structures. The same 
comments apply to the kinetic generator structure and the deformations of 
the electromagnetic fields and vectorpotentials and indeed the fluctuation 
metric itself. This completes our discussion of the parameterization of the 
hybrid fluid kinetic theory. 



3.1.1 Transformation properties of electromagnetism 

To be complete we also have to formulate Maxwell's equations in such a 
way that they are invariant with respect to diffeomorphisms in four space 
which fixes time and transform in a natural way. The standard formulation 
of Maxwell's equations in four space is with respect to a fixed metric 

0, (69) 
4ttj (1) , 

-- 3k dx k + pdt = JjW, (70) 
c 

0. (71) 

Notice that the homogenous equation transforms in a natural way with re- 
spect to diffeomorphisms in four space while the inhomogenous part of the 



* 90 d* m F {2) 



* go d *o 3 ^ 
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equation does not. We have that dF^ = di/}~ u F® = i/i^dF® = 0. We 
rectify this in a similar way as for the momentum equation by introducing 
j *ao d *go J = *gd*g = 7^ *g id d * 9id Jid- Here is the metric presented in a 
frame of reference where -det(gid) = 1 and go = ip^QM, Jid = — det(g)We 
now find the following form of inhomogenous Maxwell's equations 

* g d* g F® = 4tt]« (72) 

P(2) = F (2) /J) J(1) =J (1) /J 

* g d* g jW = 0. (73) 

When the Maxwell equations are presented in this way (or even better 
if one tensor it with the four volumeform dV g to take into account that 
it is relation between densities) they transform in the natural way under 
diffeomorphisms as 

F (2) = ^(2)^(2)^^-1*^(2)^(2)=^)^^ (74) 

T = F {2) <g> dV g = F {2) <g> dVo = F$ ® dV id , 

J = p) g, dV g = ® dV 

T = , j* = F (2) ® ^*dV ,^*dV = J(^)dV , (75) 

j« = (76) 

The upshot is that the new current and field quantities correspond to cur- 
rent density form and the field density form, J and T and they transform 
according to these.0 In case we do not have a euclidean reference met- 
ric, it can be an advantage to refer to an orthogonal frame instead with 
= F$/(-det(g ))^ and F^ = F$ /{- det(#))i It is possible to in- 
troduce electromagnetic equations with respect to J and J 7 , but it requires 
covariant derivatives which will be beyond the scope of our presentation. We 
are now ready to pullback our electromagnetic equations to the reference 
level and we find 

dFW = 0, (77) 
* go d* gQ F^ = 47rj (1) , 
*go d *g J {1) = 0. 



From a measurement point of view it also make sense to study the transformation 
properties of density forms since we always have to measure with respect to some volume 
and timeinterval. 
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The same type of equation will be fullnlled for the reference fields and cur- 
rents, F (2) and j^. If we define = ^~ U (F^ 2) + F^), F^ = F^ - 
Fq \ F^ = F^ — Fq 2 - 1 , = j'W — Jo , we find electromagnetic equations 
suitable for perturbation with respect to one diffeomorphism 

dP^ = 0, (78) 
* go d* go F^ = 47rpJ, 

In case we have multiple fluids, we would could perform the above proce- 
dure with respect to several diffeomorphisms and fluctuation metrics^ !- Notice, 
that relations like = ^~ 1 *A^ 1 \ e.t.c. are not trivial when written out 
componentwise, especially since we propose to use a fluctuating metric to 
obtain contravariant tensors. We will come back to a more elaborate study 
of the transformation and perturbation properties of electromagnetism in a 
separate paper elsewhere. 



3.1.2 Transformation properties of the stress tensor and current 
density 

Let us write our expression for the stresstensor and current density as a 
covariant twotensor and oneform respectively and study their transformation 
properties with respect to kinetic theory. 

P (2) = / ^(p (1) - mw«) ® (p« - mwW)/(«, t)d 3 p, 
J Jm 

J (1) = |^(p (1) -mw«)/(z,t)d 3 P 



W (D = u«+-A«. 

c 



We now find the pullbacked expressions 

p(2) = ^*p(2) = [ l( p (i) _ m ^(i)) g, ( p (i) - m wW)f( Zi t)d 3 p ,(79) 
J m 

j(i) = ^*j(i) = /l( p ( 1 )- m w«)/Ut)d 3 p. 
J m 

10 The same remark applies to the stresstensor which transformation properties is dis- 
cussed below. Although we will not explicitly discuss multiple fluids it is not a major 
complication from a formal point of view only for the complexity of the presentation. 
Therefore we have decided to leave this point for explicit, future applications. 
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3.2 Physical coordinates 



Our theory will not be complete before we formulate the kinetic theory in 
physical variables also. It is possible to work in canonical coordinates at 
least when there is no background, magnetic field as we have indicated in 
our work on OC theory 7 . But it requires doing the transformation from 
canonical distribution function to physical distribution function as an ac- 
tive symplectic transformation on distributions (see App.A). However, with 
a magnetic field this transformation is not so valuable since it is not a per- 
turbation. Moreover, there is the question for which coordinates one should 
specify the reference distribution. The reason why I have been reluctant to 
give up the canonical formal approach is that in physical coordinates the 
Poisson bracket itself will be perturbed. However, in section A6.3, we have 
described new physical coordinates which gives the Poisson bracket in terms 
of background magnetic fields only . We call this the interaction picture for 
the Vlasov equation. A direct extension of our work on OC 7 theory would 
be to specify the OC coordinates in interaction variables. The obvious ad- 
vantage is that now it is possible to do perturbation expansions without 
perturbing the invariant bracket and the transformation from interaction 
distribution function to physical distribution function is now really a pertur- 
bation. The relation between the interaction distribution and the physical 
and canonical distributions is given by / = (p~^*f\ f % = 4>da*f ? wnere the the 
transformation and 0^ o 1 *are described in eq.( |115| ). The relation between 
the interaction physical coordinates and canonical or physical coordinates is 
given by p = p— -Ao, p«= p — -Ax which is nothing else than the shift 
transformations described by <j)~ * and <$>~ x *. Here we use the gauge <p l = 



so Ai = —c f ~Ei(t')dt' can be given a gauge invariant meaning as propor- 
tional to the accumulated, perturbed electric field vector referred to some 
fixed time. 

The standard Vlasov equation for the distribution function, /, in eu- 
clidean physical coordinates may be rewritten in a form more suitable for 
our purposes, 



^ 7 + {/>n = o, 
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2m 

3d e 
j = /\ | B" 

dx dp p c p 
{f,H p } = J:(dpf,d P H' p ). 

Here J, {, }, A (1) and B< 2 > = dA^ are the contravariant Poissontensor in eu- 
clidean physical coordinates, the Poisson bracket in physical coordinates, the 
covariant vectorpotential, the covariant magnetic field twoform. The Pois- 
son bracket is written in standard form in eq. ([HI ). £r ,d P and d are a 
convective derivative with respect to the vectorpotential part of the electric 
field, the exterior derivative in six dimensional phase space and the exterior 
derivative in three dimensional space. In App. A we show that J =</> c *J c 
is the pushforward map of the standard canonical Poissontensor. The def- 
inition of the pushforward map is given in App. A acting on a vectorfield, 
but the action is trivially extended to higher degree contravariant tensors by 
applying the given action on each tensorindices. We have chosen to represent 
our tensors with respect to the standard basis in euclidean coordinates, i.e. 
dxi, dpi for covariant forms and -J^ , for contravariant vectorfields with 
obvious extensions for higher degree tensors. With this notation the antisym- 
metric covariant twotensor B^ 2 ) = dA^ = ^^dxj A dxk is pulled down to 
a contravariant, antisymmetric multivector identified with the corresponding 
multivector in momentum space Bp = A .[[]] Notice that with a 

timedependent vectorpotential it is not possible to treat single particle dy- 
namics as generated by a hamiltonian in six dimensional euclidean physical 
phase space except for in a convected sense. However, in eight dimensional 
extended phase space or in six dimensional canonical coordinates the dy- 
namics is described by hamiltonian generators. We show in App.B6.3 that 
the parameterization can still be presented in the new interaction picture 
by the following analogous generator equation if we take time independent 
background electromagnetic fields B ,E = — V0 O 



-j- + {f\W} = 0, 

n In a noneuclidan space we would have to be careful with the metric. All this can be 
worked out with respect to an underlying Ricmannian metric space, but we have chosen 
to postpone the description of this more general theory. 
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{f\W} = Jo :{d P f\d P W} 

2m ro ' 2m 10 ' 

7o = ~ e K 

If the background electromagnetic fields are timedependent, we can trans- 
form our Vlasov equation and canonical generator equation to interaction 
coordinates by 



J - + {f,H% = 0, {f\H% = J : (d P f\d P H% 



dr 



T o dr 
d d j e dA { 

dr' dt c dt 



(i) 

°-), 



T ^-1*T 5 A 5 J_ e "R( 2 ) 

Jo = c o J C = ^ A ^ + - B op 



T 9 A ^ - C 1 ^ 

Jc " <9x <9p _ —I j - 

Here b{, 2) = dA^, Aq"\ J and J c are the reference magnetic field twoform, 
the vectorpotential oneform, the interaction Poissontensor and the canonical 
Poissontensor. 

For those who worry about the gaugeinvariance of the above equations, 
we formulate equivalent gaugeinvariant equations by taking out the potential 
in the hamiltonian generators and add a corresponding term to the convective 
derivatives. We thereby obtain a convective derivative correponding to the 
acceleration of the electric field with obvious form of the Vlasov equation 
which is gaugeinvariantQ 



j = | + J -(-eE W ), gp ,ii 

12 In a general metric we must study covariant derivatives and parallel translation to 
give this an invariant meaning. 
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d 8 +T ( pE^ TP ~ ( P~" Al); 



To interpret the above introduction of new noncanonical coordinates we 
now give a theorem which show in which sense we can generalize our results 
in canonical coordinates to any coordinates independent if they are noniner- 
tially, convected in phase space with respect to the canonical coordinates or 
not. Since the definition of what is a hamiltonian flow or not cannot depend 
on the coordinate system, we suggest to call the flow hamiltonian if it can 
be transferred to a a hamiltonian flow by a noncanonical (not preserving the 
Poisson tensor) map as indicated in the theorem below. 

Theorem 6 The transformation of the hamiltonian generator, ip t , in canon- 
ical coordinates given by exact oneform dpip t = w c -((^„- Id)j^.) to other non- 
canonical coordinates (in general timedependent) <p(t) : z — > Z = <p(t)(z) 
is given by the action of the corresponding inverse seven dimensional map 
which fixes time (c.f. discussion in App. A) <\>~ x {Z,t) = (4>^ 1 (t)(Z),t) as 



dpip T = u-((i/>t-Id)—), 

d d d _ , .9 

■0 = Ad(4>)ip = 4> ° 4> ° 

The corresponding noncanonical Poissontensor and Poissonbracket are given 
by 

J = 4>* J c , 

{f,g}n = J : (d P f,d P g) . 

Moreover, if ift is a canonical, Poisson preserving map with respect J c , {, } 
such that ip*J c = J c and ift~ u {f, g} = {tp~ u f, xj)~ u g} , then ij) = Ad(<f>(t))iJ> 

— — -j^ 

is Poisson preserving with respect to J, {, } n such that J = J and ip {/, g) n = 

{i>~ 1 *f,i>~ 1 *g}n- 
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Proof. Let us asume that the vectorfield X^ t = J c • (d P tp t ) = (ip^ — Id)j^ is 
hamilonian generated with generator ip t . Then transform the exterior deriva- 
tive of the hamiltonian by pullback of the above map either as a result of an 
active transformation of the oneform or a passive coordinate transformation 
(c.f. the discussion below) is given by 

(j)~ u dpip t = d P ((j)~ u i) t ) = d P iJ T 

= 4- u MW. - id)^)) = ufaty* - ^)|) 

= o ^ o 0; 1 - Id)^) = u(((Ad(<j>)1>). - Id)±-) 

d - A 9 - 9 J* 

d 

X* = (0, -Id)—, 

UJ = <p UJ C = cp u c , 

Here we are treating tensors in phase space as embedded in phase space 
extended with time such that we freely interchange the action of seven di- 
mensional map fixing time and the timedependent six dimensional map 
4> on such tensors. This will also compress notation and proofs. The ac- 
tion of tensors with even a fixed time component fixes such components, but 
give extra contributions to the phase space components through terms like 
(0^ — Id)^ . Notice that for the case that <fi is not timedependent, we have 
that X t = and = , ip T = ip t while the Poissontensor, Poisson bracket 
e.t.c. is still as prescribed above. 

We now observe that vectorfield with only phase space components 

(&-/<0^ = ** t =J-(dp^r), 

is indeed hamiltonian with respect to the Poissontensor J since J-u = Id. 
The Poisson bracket is transformed by a non Poissonpreserving map as 

<t>- U {f,9} = (f>- U (J c :(d P f,d P g) = ( l>- 1 *oJ c o ( f>*:( ( j ) - 1 *d P f, ( f>- u d P g) 
= J : (dp/, d P g) = {/, g} n , f = cf)' 1 * f , g = cp' 1 * g , 
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since = cf}~ u o T o <p* for a contravariant tensor [||]. 

If t/> is a canonical Poisson preserving map such that i/)^J c = J c . Then it 
follows that a non Poisson preserving map acts as 

=> 0,J = J, =Ad(<f>)il>. 



Remark 2 .A// t/ie formal expansions we introduced earlier is valid simply 
by replacing if) =exp(X tu ) , X w = J c • dpw — > = exp(X fD ) , = J-dpiu, 
i) t = iexp(£ w )w !t , £ w = {w, •} -> r = iexp(A>)|jf , = {«;,•}„, 
# ^ # = 0-1*^ = H o 0" 1 , / - / = = VT^/o , / = ^-^/o 

§£ + {/,#} = 0- £ + = 0; 

e)t +R, H} + H, e = 0, -> 0„ T +{0„ #} n + #, g = 

i/ere £/ie notation e. zs tzsed to tofce care o/ situations where the map <fi 
also depend on the perturbation parameter e. Therefore, analogous to the 
above notation for the time generator we have that in a phase space extended 
by both time and e 

In the case that the coordinate map does not depend the parameter e, we can 
replace e by e above. This is in agreement with the above interaction picture 
philosophy where we suggest to let the transformation with respect to canon- 
ical variables and consequently the Poisson tensor be unperturbed. After the 
transformation to the interaction picture we suggest to introduce additional 
transformations to e.g. gyrokinetic, driftkinetic or oscillation center variables 
due to adiabatic or exact symmetries in the problem. These applications and 
desciption of nonlinear perturbation theory in general is outside the scope of 
this article since it requires the introduction of noneuclidean metric. 

The hybrid fluid kinetic theory can therefore now be presented in inter- 
action physical coordinates by doing the following changes in the canonical 
theory 

r = (^)- 1 7o = (^)" 1 *°#)" 1 7o, 
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il>\(.z,t) = i? lT , o (z,t)+i? 2T , n (x,t) = p • u + ^ , (80) 



■q2 

^2r(, = -7Wy--Ai-U, (81) 
lfc(z, t) = ^ le (z, ^ + ^ 2e ( X > *) = P • A + i>2e , 

VV;, = P • u , i>u = P ' • (82) 

The hybrid fluid kinetic theory can now be repeated and we find analo- 
gously to our earlier derivation in canonical coordinates 

ft* fP ,V (P-™Q 2 

r = (^y u fo, 

e . 

wi = uH A x . 

cm 

The fluid part of the hybrid fluid kinetic theory are of course the same as 
before while in interaction physical coordinates we have the resonant particle 
kinetic equation with the above parameterization given as 

dP - ~ ■ 

-L + {f,H% = 0. (84) 
ar 

The corresponding action principle for the resonant particle part will follow 
trivially from App. C by doing the above replacements in interaction physical 
coordinates. 

We are now ready to replace the above fluid symplectic and resonant 
particle generators by generators in laboratory frame which are more suitable 
for perturbations and linearization which it is our aim to do. We write 
our new parameterization as f l = (■?/> )~ u f l = (?/> ff ■ Here f l is 

the resonant particle distribution in the laboratory frame defined such that 
/ pd 3 -p =p (^ + V • (u p ) =0) . The discussion of the properties of this 
representation of the follows trivially from what we have proved before. The 
fluid kinetic generator we now define as 

= t>-^t + ^2r'M,t), (85) 
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^(x,t) = -(-^ t 2 + -A 1 -^ t ), (86) 

The resonant particle distribution in the laboratory frame is then de- 
scribed by the following Liouville equation 

§^ + {/*,ir} n0 = 0, (87) 
cLt 

-j^T + {./o> -^olnO = 0, 



& = ^+^-(^r 1 *^, 

(p-^Ai-^) 2 

The above form of the equations and parameterization are the one which is 
suitable for perturbation since everything is now really expressed in inter- 
action physical coordinates and a near identity, symplectic fluid generator 
which is directly related to the fluid generator i\) t . Moreover, there is no 
reason to change the fluid part of the theory earlier developped since it is 
already expressed in terms of the fluid generator. 

We can relate our new kinetic fluid generator to the earlier one by putting 

■0 = *0 o *0 . We observe that for the reference state t/> = 4> . We now find 
that one can express 



= (p-(mV' t + -A 1 ))-V*uo-y(^*u ) 2 , 

^ = + V'e = -^L( X > + V4( X > t) = °- 

The reason why we can take ip e = ip e is that the fluid compatibility relations 
for if> t ,ip e or ip t ,u are consistent with the kinetic compatibility relations 
^e,r' + {^e,^r' } -^ T > ,e = and i> e y a + {J e , J T& } - = 0. Also the fluid 
description of can n °t depend on which momentum coordinates we are 
using. 
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Some remarks about the form of the action principle developped in App.B 
with canonical coordinates is needed. The fluid part of the action principle, 
Af , can be kept unchanged. In interaction physical coordinates the resonant 
particle action principle in the fluid frame takes a form analogous to the one 
in App.C 

A r = J r&i, + (^r 1 *^ - ir)tfzdt = J - ^m^Mt. m 

We can develop an extended action principle for interaction physical coor- 
dinates in the way we did in App.C or equivalently in terms of the new 
laboratory frame canonical generators. The variation which gives the stress 
tensor for the fluid theory still has to be performed with respect to the fluid- 
kinetic action term Ai = — J f 1 H l d 6 zdt which we observe correspond to the 
negative time integrated internal energy. However, we can vary the resonant 
particle distribution with respect to the parameterization in the laboratory 
frame if we so wish. A variational equivalent form is therefore if one explicitly 
want to use the laboratory frame representation 

A r = J f (^ -H*- &,yzdt . (89) 

4 Conclusion 

We have described how to parameterize the solutions of the Vlasov-Maxwell 
equations by using canonical transformations with respect to a reference 
state. The canonical transformations are described by an equation for the 
timelike hamiltonian generator relating it to the local hamiltonian and the 
reference hamiltonian transformed by the Poisson preserving transforma- 
tions. Our main emphasis is on that the solutions of the Vlasov equation 
has got a composition principle for the transformations parameterizing the 
solutions with respect to a given reference solution of the Vlasov equation. In 
particular there is a specific infinite dimensional group of transformations on 
phase space extended by time- i.e. a pseudogroup, leaving the reference dis- 
tribution invariant. To see that the Vlasov equation defines a pseudogroup 
one has to choose transformations on space/time which fixes time , i.e. a 
family of canonical transformations on phase space parametrized by time. 
The pseudogroup which parameterize the solution of the Vlasov equation is 
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then given by the pseudogroup of all canonical transformations defined mod- 
ulo the transformations which leave the particle density invariant. We also 
show how Maxwell's equations can be formulated such that it transforms in 
a natural way with respect to diffeomorphisms. 

We suggest that the composition principle coming from the underlying 
pseudogroup are of fundamental importance. The principle is such that it 
is possible to specify new a priori information in the mathematical structure 
of the Vlasov equation. Although, the introduction of such a priori infor- 
mation constrains the experimental situations and the physical processes for 
which the theory is applicable, the theory has a gain in the possibility of 
modelling kinetic processes in complicated background states. The model is 
such that any perturbation theory based on it will preserve the number of 
particles in an invariant way since we are basically using mathematical en- 
tities which are not coordinate dependent. Our philosophy is therefore that 
in any modelling effort such introduction of new a priori information is a 
necessary first step. In fact, since we have replaced the continuity equation 
by a constraint equation for Poisson maps on phase space parameterized by 
time, the introduction of such a priori information will change the structure 
of our equations. Here we specifically use the composition principle to make 
the Vlasov-Maxwell equations into a hybrid fluid kinetic theory and in ad- 
dition to formulate the equations for the divergent and rotational modes in 
any background. (This application by no means exhaust the vast number 
of applications for this fundamental principle in averaging, separation and 
perturbation techniques. We are currently exploring a few of these both in 
plasma and fluid theory.) The success of the technique in this case is due to 
that we are able to identify a canonical transformation which after integra- 
tion over the momentum coordinates coincides with the parameterization of 
the fluid density. By this trick we are able to define new particle density co- 
ordinates and corresponding parameterization which contains no deviation 
on average from the density and momentum of the reference state. How- 
ever, the new density contains resonant particle effects and higher order fluid 
moments. We succeed in presenting the theory in new physical coordinates 
which we call interaction physical coordinates and the fluid kinetic generators 
in a form which is such that it is suitable for near identity transformations. 
Thus we have obtained a theory where kinetic and fluid effects are naturally 
separated. 

Such a decomposition of kinetic theory into hybrid fluid kinetic theory 
is bound to influence the way one define and think about physical phenom- 
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ena in plasmas like resonant particle distribution, dispersion relation, wave- 
resonant particle distribution effects, instability e.t.c. Work are under way 
to separate higher order fluid moments like stress. The linearized equations 
can be presented by Hermitian operators with respect to an indefinite inner 
product, but now with a mixed fluid kinetic inner product. If the system con- 
tain some extra symmetries (exact or approximate) either in the reference 
state or in an intermediate state it might be convenient to introduce other 
coordinates than the interaction physical, euclidean coordinates we have in- 
troduced. Examples of such symmetries are gyrophase and wave phase sym- 
metry or simply an ignorable coordinate in the plasma description. This is 
outside the scope of the present work since one has to think carefully about 
noneuclidean, invariant descriptions for kinetic theory. 

5 Appendix A The pseudogroup connected 
to continuity equations 

Pseudogroups are infinite dimensional groups. We will not give the detailed 
definition of infinite dimensional Lie groups since we would then have to in- 
troduce much more mathematical machinery than we intend to do here, but 
the interested reader may find it in the references (5,6). For our purpose we 
will only notice that some equations in physics define two type pseudogroups: 
one type which fixes the physical field in question and another one which de- 
forms and parameterize solutions of the the corresponding physical equation. 
Mathematically this means that the equations in question allows an infinite 
dimensional symmetry taking solutions into solutions. 

We define pseudogroups with smooth structure in the following way: 
Def. The totality T , of smooth maps on a space M form a pseudogroup 
I' {o : M — M o e C°°(M, M)} if 

i) / o g e r when /, g e V , 

ii) id , the identity map is in T , 

iii) f og- 1 eT when f,geT. 
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The continuity equation in fluid mechanics 

| + V.(v„) = 

, can be replaced by : 

Theor. The continuity equation can be parameterized by a family (with 
respect to timeparameter t) of transformations on space M , ij)(t) : M — > M 

p(t) = m*Po(t) = J(t)Po(t) ° (90) 

Here J(t) =| — | is the Jacobian of the map fp^it) ■ The velocity 
field defines an equation for the parameterization with respect to the reference 
velocity field by 



v(t) = il>(t) • v (t) + k = ij} t + ^(t).v (t) + k (91) 
V • (kp(t)) = , k = V(t)*k , (92) 

^(f).v (t) = (v • V^t)W(t)- 1 . 
Proof: 

Here the map i/)(t)* defined above is the standard pushforward map. We 
parameterized the density at a shifted time by 

p(t + s) = J^ Q (t,s)o^(t)(^o(t, s)*- 1 o V(t)* _1 )p (t + s) , 
where we have decomposed ij)(t + s) = if) (t,s) o if)(t) so that if) t = 

ds 



=0 . The following relations follows 

— ^— | s=0 - V tf t , 

J^ Q {t,s)o^{t) = (Jv,(t) o V (^s) _1 )^o(M) ' 

(9 

p(t)otj> (t,s) 1 I s=0 = -Vt • Vp(t) , 



(9s 

<9p (t + s) 



(9s 

These relations immediately gives us that 
dp(t) dp(t + s) 



s =o = -V • (vop (i)) . 



<9t ds 



o= -p(t)V ■ip t -ip t - Vp(t) - • (V • (vop (t))). 
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Since we have the identity 

^)" 1 *(V-u) = ^V.(J^( t) ^(0.u), 

-V • (vp(t)) = -V • ((^ + V(t)*v )p(t)) , 
ip t + ?/>0)* v o + k . 

Here k is any vectorfield such that V • (kp(i)) = 0. Such vectorfields can 
also be parameterized with respect to the reference density according to the 
above identity for divergences 

k = VW*k , 
V-(koPoW) = 0. 

In the following we will often not write this additional freedom explicit. 
End of proof. 

We see that in the case that the reference velocity is zero the above 
parametrization is equivalent to the Lagrangian description of fluids , but 
transported back to the Eulerian velocity v(t) by the inverse mapping. One 
can verify that the parameterization of the continuity equation is compatible 
with composition of families of smooth maps since 

J<P(t)orp(t) = (Jil>(t) ° <t>{t) 1 )J<P(t) , (93) 

(<p(t)o<4>(t))* Po (t) = il>(t) • (<p(t) • Po (t)) , 

v(t) = ((/)(t)oV(t))«v o (t) = 0(t).^(t)«vo(t)) 

= t + 0(t)*(^ t + ^(t).Vo(t)) . 

The family of transformations defined above does not conform with the 
definition of pseudogroups because of the timedependence. However, if we 
look at our family of transformations as transformations on space-time X = 
M x R which fixes time 0(x,t) = (0(t)(x),t) they can be viewed as member 
of a pseudogroup. Moreover, the velocity in this extended space is naturally 
defined as v = (v(-),l) while the family of density maps in three space is 
transcribed to the density p in fourspace which has the same value evaluated 



one obtains that 

dp(t) 
dt 

V = 
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at corresponding points and time. Now it is possible to express the above 
parameterizations in a more compressed form as 



p = • p = J o x ) , (94) 

v = 0*^0 = Oo ■ V0) O (f)- 1 , 

J = i 8£.| = |*g£| =J(t) . 



Here we have used V as the gradient operator both in three and four 
space. When we in addition identify the continuity equation in four space 
V • (vp) = as an infinitesimal Lie equation, it is clear how the density 
structure define a pseudogroup T d which leaves the density in four space 
invariant 5 ' 6 

T d = {0 G C°°(X, X) | . p = p , 0(x, t) = (0(t)(x), t)} . (95) 

The above pseudogroup is a finite Lie equation and the infinitesimal ver- 
sion of it corresponds to the continuity equation obtained by taking the in- 
finitesimal map = Id+ev-V .... in the finite Lie equation. The parametriza- 
tion of the solution space of the continuity equation as given by eq. ( |94| ) with 



respect to a reference distribution is generated by the pseudogroup of all 
smooth transformations on space-time which fix time T t = {0 G C°°(X, X) \ 
0(x,t) = (</>W(x),t)}- 

Def. We define the density leaf fixed by a density p as Po = {p I P = 
0«p o ,V0 G T t } . 

We notice that it is also possible to define a pseudogroup of smooth trans- 
formations which leaves p invariant r d which in fact can be transported to 
every element in the density leaf by T p d = {-0 = 0o-0o0 _1 | i/j e T d , G r t }. 
In this sense it is possible to look at T d as a type of generalized gaugetrans- 
formations with respect to a given leaf and moreover it will be sufficient to 
generate the leaf by f = T t modulor d with the above transportation. (In fact 
it corresponds to that one has to introduce a semidirect product as group 
product in T.) We will not pursue the generalization of gauge transforma- 
tions on the density leaf any further in this paper as it also will need more 
mathematical background than we intend to show here. 
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5.1 Geometric interpretation of the fluid generator 

In our derivation of the parameterization the new vectorfield ip t appeared as 
the derivative with respect to the near identity map tj) (t, s) . This map have 
the property that ip (t, 0) = Id . Now, we shall think of our transformations 
with respect to a specific density leaf Pq given by a reference density p . Then 
it is realized that p(t + s) = ij> (t, s) • p{t). Therefore for s = 0, i/) is the 
identity map at the point p G Po • Consequently, one can also think about 
the vectorfield ip t as a point in the space of vectorfields X p at the point p. 
At the reference density we have the reference space of vectorfields X p . It 
turns out that all vectorfields can be pulled back to the space of reference 
vectorfields. This is done by defining the related near identity map ip by 

# + «)=#)oi(i,i) . 

Again we see that -^ (t,0) = Id . This time we have that 

ipit)- 1 • p(t + s) = p (t, s) = V> (t, s) • p (t) . 

Therefore one realizes that the corresponding vectorfield i}> t = 9 ^q^' s ^ \ s =q£ 
X Pq . Moreover, it is verified that 

^ = vW*^ = (^-vv(t))oVW- 1 , (96) 

v = #M^ + v ). 

Notice that a sum of vectorfields at the reference space of vectorfields ordered 
according to the connected composition, corresponds to a total vectorfield at 
P by 



i=i 

n 
i=l 

n-1 

*l>t = *l>t,l + ^l(*)*^t,2 + - + (J! ^t(*)*°)^t,n 

i=l 

n-1 
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This is in fact our fundamental relation which tells us how to extend the above 
parameterization of the solution of the continuity equation to a composition 
of several and in principle also infinitely many transformations. 



5.2 Variational relations and perturbations 

We need to define perturbations and variations of our quantities with respect 
to generators in the space of vectorfields X p which as we have demonstrated 
can be pulled back to X p . We do that by first defining the perturbed family of 
maps with respect to one parameter e, tp(t, e) such that tp(t, 0) = ip(t). From 
this map we define the near identity maps ip (t, e, 5) = if)(t, e + 5) o if)(t, e) _1 
and ifr (t,e,5) = il)(t,e)~ l o if)(t,e + 5). One realizes as before that on the 
density leaf Pq these maps are deformations from densities p{t,e) and Po(t) 
respectively. This follows by defining the perturbed density p(t, e + 5) = 
i>(t,e + 5)»p (t) = if) (t,e,5)»p(t,e) and p {t,e,5) = ^{t,e)' 1 • p{t,e + 5) = 
ip (t, e, 5)»p (t). We can now define the generating vectorfields for the defined 
e perturbation as 

dij} (t,e,5) _ 
We = l<5=o e X p , (99) 

We = q- 5 U=o e ^tp • 

Instead of one parameterfamilies of deformations, one could define many 
parameter families of deformations or simply deformed maps if){t) of ip(t) 
without any reference to any parameters. We define the corresponding near 
identity maps y (t) = jp(t) o ^(t)^ 1 and V> W = w{t) 1 ° w(t) an d the 
deformed densities ~p(t) = V'o(^) *P(t)i Po(t) = V'o(^) *Po(0- The correspond- 
ing generating variations of these maps are then defined as the infinitesimal 
vectorfields 

5V = o ^{ty 1 e x p , (loo) 

5 W = w {t)- l o5 w {t)eX Pi) . 

We understand from the above that to formulate variational principles on 
a density leaf has certain nonconventional aspects due to that the fields in- 
volved are generated by the action of maps on densities. We give the following 
identities needed to do variations on density and the pushforward velocity. 
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These can be verified by doing infinitesimal variations of the corresponding 
parameterized fields. 



5p{t ) = -V-(^PW), (ioi) 
tf(^(f)*v ) = [5ip, ^(t)*v ] , 

-V-(^ e p(t,e)), 



e) 



de 

d(ij>(t,e)*v ) 



[if> e , if>(t, e)*v ] 



<9e 

Here the symbol [•, •] is the usual vectorfield bracket which is defined as 
[X, Y] = X-VY — Y- VX . The deformation generating vectorfields and the 
generating vectorfields in the timelike direction satisfy a certain compatibility 
condition in the space X p . 

Theor. When Si/j,ip e , tp t , k G X p , we have the compatibility conditions 

6il>, t -6tl> t -[6il>,tl> t ]+k = 0, (102) 

*l> e ,t-'>l>tK-[ll> e ,'*l>t] + k = 

V-(kp) = 0. (103) 
For Sip, ip e , ij> t e <Y Po , we have the compatibility conditions 

<ty, t -<ty t + [<ty,$ t ]+ko = 0, (104) 

fct -^t,6 +[^ e , + k = , 

V-(k p) = 0, k = V(t)*k . (105) 

Proof: The proof of the above result follows simply by writing out the 
compatibility conditions for the two equal variations 5j^p(t) \ p = p = -^Sp(t) \- p=p 
. If we use the above formulas, it is obtained that 

V-((<ty, t -6iJ> t -[6il>,iJ> t ])p(t)) = 0, 

from which the first identity follows. The second compatibility condition fol- 
lows similarly from that -J^p(t,e) = ^^p(t,e) . The last compatibility con- 
ditions in X p either follows from pulling back the above compatibility condi- 
tions in X p to X p or by studying the compatility condtions for the deformed 
density at p , i.e.5j^p (t) \p 0=Po = -§^Sp (t) \p 0=Po - In either case one finds that 

V-((6$, t -6fa + [6{i>,fa])p (t))=0. 
The derivation for the e— parameterized case is similar. 
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5.3 Parameterization of the Vlasov equation 

A special case of continuity equations are the Liouville equation and the 
Vlasov equation on the phase space of space and momentum. In this case we 
will have to deal with Hamiltonian vectorfields Xh and reference Hamiltonian 
vectorfields Xh ■ The conservation law for the particle density / on phase 
space in this case is given by the Vlasov equation (must be supplied by the 
definition of the Hamiltonian in question and the Maxwell's equations 8 ) 



^ + A' H • V/ = , or (106) 
f + {/,*} = 0. 

where {, } is the Poisson bracket. 

The volumeform in 6 dimensional phase space P with coordinates Z is 

given by the expression dV = Jyd 6 Z where Jy =| | is the Jacobian of the 
map (py 1 : Z — > z = $ V X (Z) to a standard system where dV = d 6 z . For our 
purposes we will only use standard Euclidean space with physical or canonical 
momentum coordinates which both will have volumeelement in the standard 
form. For physical and canonical coordinates with respect to Euclidean space 
we have since the vectorfields Xh preserve phase space volume that Jy = 1 
and V • (Xh) = . In more general coordinate systems, which is needed 
in gyrokinetic and oscillation center kinetic theory the conservation of phase 
space volume can be expressed by the conservation law 

f) T 

-^ + W-(X H Jy)=0. (107) 

If we combine these two equations we obtain the continuity equation in 
phase space for the quantity p = Jyf 

^XA + V ■ (X H Jyf) = . (108) 
In an analogous way as above we can now parameterize 



P (t) = m*p (t), mtc°°(p,p), 

Po(t) = Jvo(t)fo(t), 



(109) 
(110) 
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J v (t) = i>(t) • J vo (t) , 

fit) = /o^o^e wr'rMt), 

The Vlasov equation has additional structure for situations when the flow 
is described by Hamiltonian vectorfields. The maps which is generated by 
Hamiltonian flows are Poisson preserving maps and we must therefore take 
this into account. 

1. Define the space of Poisson preserving maps as T = {ip E C°°(P, P) \ 
r{f,9] = WW 9} for any f,g E C°°(P,P)} . 

We then have the following theorem: 
Theorem For ip(t) E T we have 



ll>(t)*X Ho = X^ t y-i Ho , (111) 

i/>(t + s) 

H = ^ + ^(t)- 1 *^o- (112) 
Remark: 

Notice that both the map ip Q and ip are identity maps for s = 0. There- 
fore we can regard the Hamiltonian vectorfield X^ as an element of the Lie 
algebra at the reference density structure corresponding to the pseudogroup 
r t restricted to Poisson preserving maps. It is not our purpose here to study 
this Lie algebra and its correspondence to our generalized gaugegroup (i.e. the 
Lie pseudogroup keping the density fixed) since it is best formulated with 
some more exact mathematical machinery available than we have presently 
assumed. 

Proof: 



X^(ty-i Ho , 
















OS 




dip (t,s) 




ds ^ U=0 ' 




il>(t)o^ (t,s) = i> (t,s) 
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Let the vectorfield ip(t)*X H0 act on a function on phase space / G C°°(P, M) 
. One can convince oneself that in this case one has the alternative expres- 
sion (see 0)for this vectorfield when one think of it as an operator acting as 
directional derivative, i.e. X(f) = (X ■ V)/ as is commonly done 

(m*x Ho )f = m*- 1 ° x Ho o m* ° / = m^mri, h }- 

Since the map is Poisson preserving we immediately get the result 

(m*x Ho )f = {f^{ty~ l H,} = (x m -i Ho )f . 

Consequently we have proven the first of the above results up to a possible 
Casimir generated vectorfield X c such that {C, g} = 0, Vg G C°°(P, R). For 
symplectic maps and canonical coordinates there are no Casimir for the Pois- 
son bracket while in general noncanonical coordinates there will be Casimirs. 
However, in our case this present no problem since we are only interested in 
functions restricted to the Poisson leaf generated by a reference f°. On such a 
leaf the Casimir is fixed and there will be no loss in generality to assume the 
above identity up to any function commuting with the Poisson leaf density 
/ related to a reference density f°. 

One can easily prove that ip (t,s) and ift (t,s) are Poisson maps since 
ip{t) and if)(t + s) are Poisson maps. The phase space functions ip t and if) t are 
the Hamiltonians corresponding to the Hamilonian vectorfields defined by 



d^ {t,s) _ dipit)' 1 oip(t + s ) 

ds 1 s=0 ~ Ws |s=u -^^ t 

From these definitions we derive that 

and therefore one deduce that X$ t = i/;(t)*X^ . Together with the first 
equality we then obtain that X^ t = ^Mt)*-^(i,) ■ ^ then follows that we 
can put ip t = ip(t)*~ 1 tjj t up to any function poisson commuting with the leaf 
density / corresponding to a reference density fo. End of proof. 

In canonical coordinates z = (x, p) , where p = p p + - A and p p , A are 
the physical momentum and vectorpotential, the Hamiltonian vectorfield is 
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given by X H = J c • d P H = (|^, — ^) in canonical, euclidean coordinates. 
The transformed Hamiltonian is given in the same form since the symplec- 
tic and the Poisson tensor does not change by canonical transformations, 

X«t)-HHo) = Jc-d P (m- U H ) = • In physical 

coordinates based on Euclidean space the symplectic tensor depends on the 
magnetic field which change also has to be specified, i.e. B — > B .Therefore 
we have that 



X T J H -( 9H ° dH ° B ° x dH °) 
^H — Jo ' "P-no — {-^ — , 7{ x — ) i 

op ax c op 

X m -i. Ho = J-dp(if>(t)- u H ) 

dj}{t)- u H d^(t)- u H B x &ijj(t)- u H 

dp <9x c <9p 

We have not given explicitly how the magnetic field changes under the 
action of the pseudogroup of smooth transformations on space time here. 
The answer to this question follows from the same infinite dimensional sym- 
metry for electromagnetic fields which are responsible for the usual gauge 
parameterizations. We will explore this more general parameterization of 
the electromagnetic fields in a forthcoming paper. 

The connection between the canonical distribution function and the phys- 
ical distribution function in physical coordinates (x, p = mv) is given by a 
canonical tranformation in the radiation gauge as 

/ = C _1 7, (H3) 
(f)' 1 * = exp(X c ) , 

c op c 

Here is the vectorfield A lifted to a oneform. Notice that the vectorfield 
X c is phase space volume and Poisson preserving, but it is not generated by 
a hamiltonian. The explicit form of this vectorfield in other cordinates will 
follow from the tranformation properties of the Poisson tensor, J. The above 
transformation is nothing else than the shift transformation from physical p p 
to canonical p. 

• Lemma The canonical transformation <p c transforms the bracket be- 
tween two functions /, g on canonical phase space to the physcical 
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bracket between the corresponding functions f,g on physical phase 
space. 

<fc U {fi9} = {f,g}n = J : (dpf,dpg). 
Here the physical bracket is given in it's standard euclidean form 



r t -i 9f dg df dg eB .df dg . 
{f,9}n = ^- ■ jt- ~ • + — ■ x . 114 
<7x p opp opp ax p c ap p ap p 

Proof. The Poisson tensor in canonical coordinates can be expressed by the 
multivector J c = ^ A . The action of the pullback map gives (?} = 

<fc u {3 c ■ (d P f,d P g)) = o J c o 0*) : (^dpf^^dpg). The exterior 

derivative operator commutes with the pullback operator <f)~ u dpf = dpf 
and for contravariant tensors (f)JT =(j)*~ 1 o Tocjf . Therefore one have that 
fc U {f, 9 } = J : (dpf, d P §), J = <P C J C = ^ A JL + - ||)^- A 

This expression is identical to the standard particle Poissontensor in 

euclidean physical phase space variables given above. 
■. 

We immidiately notice two major problems with this bracket. It is not 
compatible with the reference distribution, f° since the Vlasov equation for 
that has to be expressed with respect to background electromagnetic fields. 
Secondly, it is not compatible with perturbation theory either since then 
one would have to do a perturbation expansion of the bracket itself. This 
completely destroys the ideas we advocated for above using canonical fixed 
brackets as a tool for invariant expansions. To resolve this in our opinion 
fundamental problem in plasma physics, we suggest to define a new physical 
distribution function given by the background fields f l = <Pco*f ■ The bracket 
for these kind of distribution functions are now transformed to the same 
form as in eq.( |114j ), but with B — > Bq. This distribution function is still a 



gaugeinvariant distribution function since we will use the gauge, <p l = where 
t 

Ai = — c J Eii(t')dt' has a physical meaning in terms of the timeintegrated 
perturbed electric field, Ei = E — Eo . The relation between the physical 
distribution function and the interaction distribution function is given by 



/ = <t>;i*r, (us) 
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(Pel*' = eX P(^ Cl ), 

X cl = -Ax • A = -Jo • (-AW)) (116) 
c ap c 

Jo is the euclidean physical coordinates Poisson tensor with B — > B . It 
seems fitting to call this description of the Vlasov fields the interaction picture 
since it is now possible to separate background and fluctuating quantities in 
an invariant way suitable for perturbation theory. 



6 Appendix B Hybrid fluid-kinetic action prin- 
ciple 

We will in this appendix study the action principle for the hybrid fluid-kinetic 
theory. In two other works 4,8 we have elaborated on the action principles for 
the Vlasov equation and the ideal fluid equations respectively. Our approach 
is based on varying the generators of the underlying infinite dimensional 
group acting on the respective densities. The basic method is quite different 
from the approach of Larsson 1 ' 2 which is using canonical conjugate variables 
on the accessible leaf. However, our method can be revised to introduce 
canonical conjugate variables with certain differences since our action also 
explicitly takes into account the group composition law and the compatibility 
conditions. The action principle for the Vlasov equation is (the Maxwell 
equation has it's own action principle which in fact also can be parameterized 
by an infinite dimensional group) 

A p = J f(ip t - H)d 6 zdt. (117) 

The compatibility condition for the parameterized phase space density can 
then be used to formulate a revised action principle for densities which de- 
pend on an additional formal perturbation parameter e, i.e. f(t,e) 

Ai 1] = J J f& e ,t + R, H} - HJdPzdtde (118) 

In eq. |117| we treat / as parameterized by symplectic transformations with 
respect to a reference state f°. The variation is nonstandard in the sence 
that the action is varied and sought stationary with respect to the the in- 
finitesimal Hamiltonian generator Sip such that 5f = {Sip, /}. (However, this 
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action imply standard variational principles by the introduction of the re- 
vised variational principle through one parameter variations.) Moreover, the 
variation of the generator ip t is determined through the compatibility relation 
for the variation 5^ = j^Sf which leads to the compatibility condition for 
variations 



5ip, t -Sip t + {Sip,ip t } = OmodA;, (119) 
k = if)~ u k°, {k,f} = 0. 

The function k has no influence on the variations. Note that if the variations 
are restricted to a one parameter group, the above compatibility condition 
is equivalent to the one we have derived before since then Sip = ip e 5e . This 
means that for both action principles one obtain the Vlasov equation by the 
variations 

In the revised action principle it is possible to introduce a variation with 
respect to 5f keeping ip e fixed since variation with respect to the one pa- 
rameter generator ip e is only a subvariation. The variation with respect to / 
then gives the compatibility condition as one of the variational equations. 

5A (1) 

= V M - H, t +R, H} = . 

In fact, there is no reason why one could not introduce many (even infinite) 
parameter groups if this is suitable for the problem at hand. If we want, 
we could also give up the explicit parameterization of / through symplectic 
transformations in the revised action principle and formulate a canonical 
field theory for canonical conjugate variables ip e , f as Larsson 1 ' 2 do. In this 
case one would have to introduce f, £ = {ip e , /} as an additional constraint. 
In our paper 4 we do this by introducing a Lagrange multiplier, but also by 
embedding the problem in a larger double symplectic space which includes 
also the e -dynamics. For some purposes this might be a somewhat restrictive 
point of view, e.g. if one want to derive model equations based on several 
layers of transformations as we want to do. 

We are now in a position to formulate a new hybrid fluid kinetic action 
principle where we restrict one part of the symplectomorphism, tp, to corre- 
spond to what we have found in section 3 to be equivalent to volume density 
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preserving transformations in space (both parameterized by time even if we 
do not explicitly indicate it). The second part of the composition corresponds 
to an incoherent kinetic transformation ip due to higher order Hamiltonian 
generators than linear in the momentum coordinate. The total hybrid fluid- 
kinetic action restricted to such a composition takes the form 



A H = Ar + A F , 



A r = / f(ijj t + if, rj) t - H)(Fzdt = / V* - H)d b zdt 



Here the Hamiltonian H and the related H is defined in eq.( 29). Up 
to variational equivalence (which after all is what is important in a vari- 
ational principle), we can freely move the action of a symplectic transfor- 
mation between a phase space density / = <j> f with suitable decaying 
properties in infinity and a multiplying phase space function J fgd 6 zdt -<=>- 

J f((j)*g)d 6 zdt. Therefore equivalently the incoherent part of the action prin- 
ciple restricted to fluid orbits can be written 

A r = J f(jp t -^*H)d 6 zdt (120) 

Such changes between variational equivalent forms of variational principles 
will later on be freely done without further mentioning. The fluid part of the 
action principle has the form 8 (here the density parameterization is given in 
the action principle) 

A F = j p(wfi-y)d 4 x, (121) 

e . 

w = u H A , 

mc 

P = VV, (122) 
u = V't + ^uo. (123) 



Remark 3 The fluid action, A F = Af + Ai, can further devided into one 
part which is simply the momentum space integrated Af — J (J fiptd 3 p)d A x 

-—2 2 

if we identify \ with ^ and one part which could be identified as the elec- 
tromagnetic/fluid interaction part, Ai = j • ud^x. 
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Lemma 7 The above fluid action is again nonstandard in the sense that the 
variations has to be done with respect to a infinitesimal variation of the fluid 
generator Sip for the quantities which are parameterized with respect to the 
reference fluid state p°, u° 

5 P = -V-(Sipp), (124) 
Sip, t = Sip t + [Sip, ip t ] , 
Su = Sip, t -[5ip,u] . 

The variations with respect to u and the electromagnetic potential are stan- 
dard. We have in our earlier work 8 , found by using the above relations that 

SA F ,du 1 

dip at m 

f L = — (E+— u x B) , 

m c 

E = B =V x A , 

c at 



5A 



F 



5u 



p(u — u) = 0^>u = u. 



It also seems natural to call the term Ai = — J fHd % zdt which correspond 
to the internal energy for the fluid-kinetic interaction part of the action. The 
variation of this part of the action with respect to the fluid generator gives the 
divergence of the stresstensor needed to complete the fluid momentum equa- 
tion. In this formulation of the action principle the mass density continuity 
equation is implicitly given by parameterization of density. 

5 A! f 5ip~ 



Sip 



J ^{f, H}d 3 p = -J p{f, H}d 3 p = -V • (P) . 



Here we have used the obvious lemma valid for phase space densities with a 
suitable decay in infinity and an appropriate class of phase space observables 
which g(z,t) belongs to 



Lemma 8 / G(x, t){g, /}rf 3 p = - V • ( / g/d 3 p)G(x, t). 



54 



The variation with respect to the infinitesimal generator 5ip gives us the 
reduced Liouville equation on fluid orbits with respect to similar variational 
rules as we discussed in eq.( |119| ) 

^ = _f_ (7,^ = 0. (125 , 

6.1 Revised hybrid fluid-kinetic action principle 

From the parameterized version of the hybrid fluid-kinetic action principle 
it is possible to derive a revised action principle in the same way as we did 
above for the Vlasov action principle. This is done simply by assuming that 
/, /, ij> and ip depend on an additional formal parameter e. We then find 
the revised hybrid fluid-kinetic action principle 



1 

AS = J |(^-(-p(^ + fi-Vu-if x -V 1 P) (126) 
o 

-(^•w)(^ + V-(up)) (127) 

+p(u, £ + V £ • Vu) • (u - u)+p— A,, • u)d 4 xde (128) 

cm 

i ^ 
+ // 7fc + ft + TH} - ^^(Fzdtde . 
o 

The underlined terms in the fluid and kinetic part of the action are inter- 
changeable forms of the same term. With this revised action principle we 
obtain the same equations as above by varying with respect to t/> e , u, e and 



7 App. C Rotation and divergence defined 
in an invariant way 

One way to define the rotation and divergence of a fluid element in an invari- 
ant way is through the Hodge star operation relative to a metric [||]. Another 
more intrinsic way is through the Lie derivative of a volume element. Our 
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interest in this is motivated by the the need to formulate physical equations 
and here fluid dynamics in such a way that they transform naturally with 
respect to diffeomorphisms. From a practical point of view this is needed to 
formulate perturbation theory with low complexity and new models. How- 
ever, from a more fundamental point of view there is a need for an intrinsic 
description of observable quantities like rotation and divergence of the flow 
of a fluid element. The definition of the Hodge star operator with respect to 
an invariant volume element dV = JdVo = ip~ *dVo, give that one can easily 
check that it must behave naturally with respect to diffeomorphisms 

aA* s /3 = g-\a,p)dV =< a,g _1 (/3) > dV (129) 
* g = f^o^of, (130) 

g = V _1 *go- 

Here * g and * g0 is the Hodge star operation with respect to g and g respec- 
tively and a and f3 are forms in A k T*M for some k (k — 1, 2, 3 for us in three 
dimensional space). By abuse of notation we will use the same notation for 
the metric g and its inverse g _1 as for the maps induced by them, e.g. here 
g _1 /3 = g _1 (-,/3) is a contravariant multivectorfield in A k T M. Moreover, 
< •,• > is the standard contraction between A k T*M and it's dual space 
A k TM. A more direct description of the action of the Hodge star operation 
is formulated by 

Lemma 

*s v = Jg ° So 1 ° *«>*>, 

where v is a form in A k T*M ,k = 1, 2, 3. 

Proof. With respect to the metric go we have that a A * g0( 5 = 
< a,g 1 (f3) > dV . Therefore one find that a A * s v =< a,g " 1 ((g o 
g-» > dV =< a, g \J(g o g-» > dV 

= a A * go (^(go ° g _1 )f ). Since this equality is valid for all k-forms a, the 
first part of the lemma is proved. For the second part of the lemma we use 
that we can write v = * s v' = * so v' where v' = * s v and v' = * so v. Here 
we have used that in three space (same equality up to sign in some other 
dimension) * g o * g = * g0 o * g0 = 1. The first part of the lemma then implies 
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that v' Q = * go v = J(g o g _1 )f' and consequently v' = * g v = j(g o gg x ) * g0 v 
. This proves the second part of the lemma. 

We are now in a position to state Hodge decomposition (we do not con- 
sider singular contributions) with respect to a general metric as 

Theorem 9 Forv, A, rj as forms in A k T*M, A 1 'I' M. K n - k ~ l T*M respec- 
tively (n = 3 for M three dimensional) , we have that 

I 

v = dr] + * g dA = dr] + —g o g" 1 o * So A, (131) 

Jg-\v) = Jg-^d^ + g^i^dA), 

v = ip- 1 *(df) + * So dA),ri = ip- u ri, A = -4)~ U A. 

The proof follows from direct use of the above lemma. 
We will now specialize to oneforms and define divergence and rotation 
with respect to a transformed metric g. 

Definition 10 div g (v) = * g d* g \^\ curi g {y) = g ^gcM 1 )), = g(v). 
Here v is a vectorfield in TM, v :.\I — T.\I. 

This definition leads to the following theorem 

Theorem 11 

div s (v) = jdiv So (Jv) = ij)' 1 * (div g0 (v)) , v = t/^v, (132) 

curlgV = jcurl so (g 1 o g(v)) = if) ^(curl go v) . 

Proof. We prove this by applying the above definition for divergence and 
curl. div e (y) = * g d* g v^ = * g d* g0 (Jv 1) ) = j* go rf*g (JvJ 1) ) = Jdiu g0 (Jv), vj, 1} = 

g (v). On the other hand, we have that * g d * g = ip~ u (* go d * g0 Vq 1 - 1 ) = 

if)" 1 * (div g0 (v)) , = go(v). Similarly for curl we have that curl g (v) = g _1 (*grfv^^) = 

lg-1 o go go 1 (*go d (go(go 1 ° g( v )))) 

= Jgo 1 (*m d (So(So 1 ° g( v )))) = i curl so(So 1 ° g(v)). On the other hand 
we have that g~ 1 (* g c?v( 1 )) = ?/>*(go 1 (*go^^ 1 ^)) = V'*( CMr ^ g o( v ))) an d the 
theorem is proved. 

A more geometric way to study divergence independent of metric is 
by Lie derivative of the invariant volumeform with respect to the veloc- 
ity field C{y)dV = di v dV = div g {y)dV = Jdiv s (v)dVo = di v (JdV ) = 
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di(jv)dVo = C(Jv)dVa = div so (Jv)dVo. On the other hand we have that 
C{y)dV = ^ _1 *(£(v)dVo) = ip' 1 * (div g0 (v)) JdV . We therefore observe that 
a divergence free vectorfield or purely rotational is simply a vectorfield v c 
in the kernel of the Lie operator, i.e. (by the way an infinitesimal Lie equa- 
tion by App. A) C(v c )dV = div s (v c )dV = div So (Jv c )dV = 0. This is 
consistent with a parameterization of a rotational vectorfield for an invari- 
ant volume element as v c = cur/ g (A) =jCur/ g0 A = i/^cur/g^A)) where 
A = g^ 1 o g(A), A = ifj^A. The Hodge decomposition then give us that 
the velocity field can be decomposed with respect to an invariant volume 
element as 

v« = -drj + * g dA« = -drj + -g o go 1 (*g dAW). (133) 

u 

However, with respect to the reference metric g we could consider the 
decomposition of Jv^ or Jv. In fact with respect to a reference state where 
the fluid is fixed and homogenous (pg is constant and is the Jacobian with 
respect to the corresponding diffeomorphism) , we could just as well multiply 
by the constant mass density and obtain since p = -P p§ a decomposition of 
pvW or pv as pv = — pV s i] + curl So (p{ ) gQ 1 (A^)) , V g r] = goog -1 ^). This 
decomposition is interesting since it is exactly the one we need in connection 
with the discussion of the pseudogroup defined by the continuity equation 
defined in App. A. 

7.1 Rotational bracket structure 

We will represent a rotational vectorfield v c (or one formv^^by a Pfaff de- 
composition of the one form which in the nonsingular case is A* 1 ) = 
ad(5 + d'y. The rotational vectorfield is then v?' = * s d(adf3), i.e. we can 
mod out 7. We want to think about (3 as a family of level surfaces (foliations) 
which the rotational vectorfield is situated on. With this interpretation in 
mind we will use the notation v c = = g _1 (vc 1 ' ) ) , v^ 1 ' = X« for a given 
metric g. Analogous with the Poisson bracket in phase space we define the 
new rotational bracket on the foliations defined by (3 as {a, f}/3 = — X^(/) = 
* g (X^ (1) A* g df) =< df, g _1 (x£ (1) ) > = g-\df, X^ (1) ). The equivalence of the 
first two definitions comes through that both of them are by trivial use of the 
definitions equivalent to the third expression. We can now find the follow- 
ing lemma valid for rotational one forms represented with respect to a given 
foliation. 
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Lemma 12 



*„(XC AX?) = {c^aaVfr (134) 



rv/3 v/3 1 _ v/3 

X^(/(/5)) = where f is differentiable. 

Proof. For the proof of the first identity we find that * g (X^ 1) A X„f ) = 

*g(*g(dai Ad/?) AX2 1} ) = * g (d / 3A* g (Xf 2 1) Aden)). Now, we use the identity 
* g (a« A * g (b( 1 ) A C W)) = g-^aW, c«) b« - g" 1 ^ 1 ), c«) a« (analogue to 
the classical triple crossproduct) to derive * g (X^ (1) AXq^) = — g _1 (Xa 2 1 \ dai)d(3+ 
g- 1 (X^\dp)da 1 . But we have that Xg 2 (/3) = - * g (* g (da 2 A d/3) A d(3) = 
which implies * g (X^ 1 1 ' ) A xff) = {«!, a 2 }/3 This statement also proves 
the third part of the lemma since X^(/(/?)) = /'(/3)X^(/3) = 0. The sec- 
ond identity is established by that for purely rotational vectorfields (c.f. the 
section about parametrization of hybrid fluid kinetic theory in the main text) 

PC XJJW = - * g d(* g (Xff) A X^)) = - * g d({ 0l , « 2 }d/3) = X^ i a2} . 

Therefore there is a Lie antihomomorphism between the Lie algebra of rota- 
tional vectorfields and the corresponding Lie algebra with respect to the 
rotational bracket on the space of rotational potentials for a given fam- 
ily of foliations of space. All the usual relations for Lie algebra's like Ja- 
cobi identity e.t.c. follows for the rotational bracket structure through this 
antihomorphism. ■ 
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